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1.0  INTRODUCTION 


1.1  GENERAL 

Two  broad  classes  of  continuum  theories  xppeax  to  exist  with  r^ard  to  the  modding  of  the 
loading  response  of  granular  materials.  These  are  the  phenomenological  (macr(»copic)  and 
microstructural  (particulate)  theories.  The  i^ienomenological  theories  attempt  to  devise 
ccmstitutive  relatimiships  based  on  analysis  of  iht  experimentally  determined  stress-strain 
response  of  granular  materials.  The  microstructural  theories  are  founded  upon  the  premise 
that  actual  measurement,  geometric  modeling,  or  statistical  modeling  of  the  granule  (xmtacts 
or  granular  motions  will  provide  insight  into  more  accurate  nuxiding  of  the  material 
response.  The  microstructural  theory  can  be  directly  linked  to  the  continuum  theory  if  one 
maintains  that  the  specimen  boundary  distances  (or  loaded  volume)  are  large  in  comparison 
to  the  size  of  the  granules  or  the  distance  between  them.  This  stipulation  permits  one  to 
transform  the  body  of  discrete  particles  into  a  statistically  equivalent  nuu:roscopic  !9)ecimen 
which  is  then  suitable  for  mathematical  analysis. 

Several  distinct  types  of  models  exist  within  the  two  broad  classifications  described.  The 
dominant  modeling  techniques  used  in  the  two  classy  will  be  presented  to  describe  each 
tq^noach  as  it  has  progressed  to  its  current  state  of  achievemoit.  This  will  permit  an 
objective  assessment  of  the  advantages  and  limitations  of  each  modeling  technique  and 
provide  an  understanding  of  the  contribution  of  this  research  effort.  Macroscopic  cmitinuum 
theory  will  be  presrated  first  in  deference  to  the  historical  contributions  made  by  previous 
researchers.  Microstructural  theory  will  follow  because  of  the  relatively  recent  interest  in 
this  area  of  research.  Prior  to  discussion  of  these  classes,  a  brief  qualitative  summary  of  the 
loading  response  of  soil  will  be  given. 

1.2  SOIL  MECHANICAL  BEHAVIOR 

Some  of  the  earliest  models  devel(q)ed  to  predict  the  loading  response  of  soils  were 
predicated  on  linear  elasticity  (Hooke’s  law)  and  nonlinear  elasticity.  This  was  principally 
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due  to  a  lack  of  general  understanding  of  soil  mechanical  response  but  also  due  to  the 
inal^ty  of  performing  the  large  number  of  numerical  calculations  associated  with  more 
advanced  models.  This  limitation  precluded  accurate  modeling  of  soil  loading  response, 
even  at  relatively  low  stress  levels,  because  soils  have  been  shown  to  undergo  plastic  strains 
at  very  low  levds  of  total  macroscopic  strain  (e  «  ID^).  The  criticality  of  modeling  soil 
fda^c  resptmse  has  been  given  considerable  emphasis  in  recent  models  and  some  (Ref.  1) 
have  entirely  eliminated  the  regicm  of  elastic  deformation.  In  light  of  these  develqrments,  it 
is  necessary  to  be  familiar  with  the  essential  features  of  the  mechanical  reqxmse  of  soils 
subjected  to  hydrostatic  and/or  deviatoric  stress  states.  Furthermore,  the  existence  of  pla^c 
strain  at  such  low  levels  of  total  stress  or  strain  requires  one  to  be  familiar  with  die 
fundamoital  equations  of  a  generic  work-hardening  elastic-plastic  (or  dastic-vicroplastic) 
constitutive  modd.  Because  such  models  require  the  use  of  a  flow  rule,  the  ramifications  of 
using  an  associated  or  nonassodated  flow  rule  nnodeling  the  response  of  granular 
materials  will  also  be  discussed. 

More  than  most  other  oigineering  materials,  the  constitutive  response  of  soils  subjected  to 
external  loading  is  greatly  affected  by  such  factors  as  physical  structure,  void  ratio,  d^ree 
of  void  saturation,  drainage  conditions,  current  stress  state,  loading  history,  and  loading  rate 
(Ref.  2).  Consider  the  example  of  soil  subject  to  a  hydrostatic  stress  state  (Fig.  1).  The 
pressure-volume  relation  is  nonlinear.  The  soil  undmgoes  permanent  volume  change  as 
shown  by  the  load-unload  cycle  of  deformation.  In  this  example,  the  permanoit  volumetric 
strain  is  dqrendent  on  the  initial  state  of  compaction  of  the  material,  the  degree  of  void 
saturation,  and  drainage  conditions.  Regarding  the  initial  state  of  compaction,  a  dense  soil 
will  exhibit  relatively  little  volumetric  strain.  On  the  contrary,  a  highly  flocculated  soil  can 
render  extreme  volumetric  strains.  In  addition,  soil  typically  exhibits  a  small  amount  of 
permanent  volumetric  strain  during  the  unload-reload  cycle. 

The  shearing  stress-strain  response  of  soils  is  much  more  complex  than  that  of  simple 
hydrostatic  response.  Figure  2  shows  typical  qualitative  responses  of  sheared  soils.  Hiis 
figure  indicates  the  influence  of  confining  pressure  (hydrostatic  stress,  aj  on  the  shear 
strength  of  the  soil.  Goierally,  increased  confining  pressure  results  in  increased  stwar 
strength. 
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Figure  1.  Typical  behavior  of  soil  under  isotrqnc  loading 
conditims  (from  Ref.  3). 


Figure  2.  Typical  triaxial  stress-strain  curves  for  soil  (from  Ref.  3). 
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Soils  ty{»caUy  exhibit  volumetric  changes  when  subject  to  a  shearing  (deviatoric)  state  of 
stress  or  deformation  (Fig.  3a).  The  volume  change  can  be  either  cmitractive  or  dilative 
depending  on  the  initial  state  of  compaction.  A  loose  sand  or  normally  ctmsolidated  clay 
will  usually  contract  under  shear  whereas  a  dmse  sand  or  overconsolidated  clay  will  tend  to 
dilate.  Figure  3b  shows  the  responses  of  strain  softening  and  strain  hardening  which  can  be 
associated  vdth  dilative  materials.  A  strain  softening  (dilative)  material  shows  a  relatively 
well-defined  peak  stress  at  relatively  low  strain.  Continued  increase  in  implied  strain  results 
in  a  decrease  in  the  strength  of  the  material.  This  is  a  brittle-type  failure.  Altemadively,  a 
strain  hardening  (contractive)  material  will  undergo  an  increase  in  stress  (strmtgth)  with 
strain  until  some  peak  stress  (strength)  is  achieved.  Figure  3c  provides  the  same  informaticm 
as  Figure  3b,  but  the  response  is  expressed  in  terms  of  void  ratio. 

The  stress-strain  response  for  soils  has  been  shown  to  be  dependmt  on  the  applied  confining 
stress.  Likewise,  there  is  a  dependence  on  time  rate  of  loading.  However,  the  dependence 
is  not  so  clear  cut  and  varies  with  respect  to  the  frictional  and  cohesive  properties  of  the 
material.  To  clarify  this  point,  the  limiting  cases  of  purely  frictional  and  purely  cohesive 
materials  will  be  presented.  Generally,  the  shear  strength  of  frictional  materials  is  primarily 
dependent  on  the  confining  pressure  and  relativdy  independent  of  the  strain  (or  shearing) 
rate.  Alternatively,  the  shear  stroigth  of  a  a)hesive  material  is  primarily  dependoit  on 
strain  rate  and  much  less  so  on  confining  pressure.  The  shear  strength  of  mixed  soils  (those 
having  both  frictional  and  cohesive  attiibutes)  depends  on  both  the  confining  pressure  and 
shear  rate.  The  effect  of  shear  rate  on  mobilized  soil  strength  is  shown  qualitatively  in 
Figure  4,  which  depicts  the  results  of  four  axisymmetric  triaxial  tests.  In  this  test  smies, 
four  similar  specimens  have  been  isotropically  consolidated  to  the  same  confining  pressure. 
One  of  the  specimens  is  then  sheared  by  static  load  application  and  the  other  three  are 
sheared  at  different  rates.  The  key  feature  apparent  in  this  figure  is  the  direct  relationship 
which  exists  between  shear  rate  and  material  strraigth.  The  strength  of  the  material  increases 
as  the  shear  rate  increases.  The  strength,  therefore,  is  defined  by  an  envelope  known  as  the 
"dynamic  failure  envelope.”  This  soil  attribute  permits  a  specimen  to  exhibit  strain 
hardiniing  under  one  loading  condition  and  strain  softening  under  another.  Figure  S  shows 
the  results  of  axisymmetric  triaxial  tests  being  conducted  on  two  similar  specimens. 
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O  dense  sand  or  OVERCONSOLIDATED  CLAY 


(2)  LOOSE  SAND  OR  NORMALLY  CONSOLIDATED  CLAY 


0  DENSE  SAND  OR  OVERCONSOLIDATED  CLAY 


0  LOOSE  SAND  OR  NORMALLY  CONSOLIDATED  CLAY 


0  LOOSE  SAND  OR  NORMALLY  CONSOLIDATED  CLAY 


0  DENSE  SAND  OR  OVERCONSOLIDATED  CLAY 


SHEAR  STRAIN 


Figure  3.  Typical  behavior  of  soils  tested  under  drained  triaxial  test 
conditions  (from  Ref.  3). 
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MEAN  NORMAL  STRESS,  P 


(a) 


(b) 


Figure  4.  Typical  stress-strain  curves  for  different  strain  rates 
from  an  axisymmetric  triaxial  test  (from  Ref.  3). 


O  axial  STRAIN, 


Figure  5.  Dynamic  stress-strmn  curve  for  a  variable  strain 

rate  from  an  axisymmetric  triaxial  test  (from  Ref.  3). 
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The  bottom  curve  shows  the  results  of  a  test  in  which  the  specimen  was  statically  loaded. 
The  upper  curve  shows  the  results  of  a  test  in  which  the  specimoi  was  r^idly  sheared  to 
some  peak  stress  and  then  slowly  sheared  to  some  ultimate  Oower)  diear  stress.  The 
differing  respmises  are  induced  solely  by  varying  the  shear  rate  during  dynamic  loading. 

All  of  the  information  presented  thus  far  reflects  tests  conducted  on  specimofis  in  which  the 
pore  fluid  (if  any)  has  been  able  to  drain  during  testing.  Tests  conducted  on  saturated  soils 
under  undrained  conditions  provide  results  which  differ  considerably  from  the  results  of 
similar  tests  on  drained  specimens.  The  limiting  cases  in  undrained  tests  are  those 
conducted  either  on  saturated  overconsolidated  clays  or  very  dense  sands  and  tests  on 
saturated  normally  consolidated  clays  or  very  loose  sands.  Figure  6  shows  the  qualitative 
results  of  tests  c<mducted  on  "limiting  case"  specimens  (designated  2-3  for  a  very  loose  sand 
or  normally  consolidated  clay  and  2-S  for  a  very  dense  sand  or  ovmconsolidated  clay)  in 
undiained  shear  conditions.  The  specimen  designated  2-4  has  material  properties  somewhere 
between  these  limiting  cases.  Each  of  the  specimens  was  isotropically  consolidated  to  the 
same  ctmfining  stress  and  then  sheared  (undrained)  at  the  same  constant  rate.  The  radial 
stress,  remained  constant  throughout  the  test.  Furthermore,  each  specimen  was  subjected 
to  a  loading  r^ime  which  followed  the  same  total  stress  path.  Each  specimen  experienced  a 
different  effective  stress  path. 

Effective  stress  in  a  fully  saturated  soil  element  is  defined  as  the  total  stress  minus  the  pore 
pressure. 


*  o  -  U  or  o  =  +  U 

where  a'  equals  effective  stress,  a  equals  total  stress,  and  U  equals  pore  pressure. 
Comparing  Figures  6a  and  6b,  one  can  see  that  the  effective  stress  is  the  only  portion  of  the 
total  stress  which  affects  soil  shear  strength.  Increasing  the  effective  stress  causes  soil 
particles  to  shift  into  a  draser  pacldng,  whereas  equal  increases  in  the  total  stress  and  pore 
pressure  leave  the  effective  stress  unchanged.  In  the  latter  case,  there  is  little  or  no  effect 
on  the  particle  pacldng  (Refs.  2,  4,  and  5). 


7 


VOID  RATIO  ^  SHEAR  STRESS 


PL-TR--92-1064 


PL-TR~92-1064 

Previous  qualitative  results  have  been  predicated  cm  the  results  of  tests  conducted  on 
isotrt^  sc^.  Aniso.  opy,  however,  mu  *  be  ccmsidered  as  it  is  a  ccmimon  soil  property. 
Figure  7  shows  qualitative  results  associated  with  the  conduct  of  hydrostatic  loading  tests  in 
an  isotropic  qjedmen  ami  an  anisotropic  specimen.  The  isotropic  ^)ecimen  exhitnts  equal 
strains  in  the  vertical  and  horizcmtal  directicms  (i.e.,  e, .  e,).  The  anisotropic  qwcimen 
exhibits  differait  strains  in  the  vertical  and  horuxmtal  directimis. 

The  specialized  case  of  uniaxial  strain  (e^  is  constrained  to  be  0)  must  be  considered  because 
many  tyjncal  loading  regimes  involve  the  a^lication  of  uniaxial  stress  or  strain.  Figure  8 
shows  the  qualitative  results  of  the  conduct  of  a  uniaxial  strain  test  cm  a  drained  ^lecimen 
with  typical  behavior  (Fig.  8a)  and  the  associated  stress  path  (Fig.  8b).  Figure  8a  dqncts 
the  typical  relatively  high  initial  loading  modulus  due  to  natural  cementation,  but  the 
material  subsequently  densifies  and  the  response  approaches  that  obsoved  under  drained 
hydrostatic  loading.  Upon  unloading,  the  material  exhibits  a  rdativdy  stiff  reqxmse  due  to 
interparticle  locking;  however,  this  resptmse  gradually  decreases  as  unloading  proceeds. 
Figure  8b  indicates  the  presence  of  locked-in  horizontal  stresses  ailer  the  specimen  has  been 
unloaded. 

1.3  PHENOMENOLOGICAL  (MACRQSCOPIO  CONTINUUM  THEORY 

Accurate  prediction  of  the  soil  loading  response  requires  two  principal  ingredients.  One 
must  have  a  modd  that  exhibits  the  major  features  of  the  material  response  under  three- 
dimek.  .’anal  loading  as  well  as  quantitative  determinatimis  of  the  material  parameters  which 
characterize  the  model.  It  is  critically  important  to  view  the  chosen  modd  and  determine 
the  required  matnial  parameters  as  codependoit  aspects  of  the  problem  definition.  This  will 
preclude  the  use  of  an  overly  complex  model  or  determination  of  parameters  which  have  no 
direct  influoice  on  the  accuracy  of  the  prediction.  It  is  implied  that  only  those  testing 
devices  ^ch  subject  a  spedmen  to  arbitrary  three-dimensional  loadings  (Ref.  6)  are 
suitable  for  obtaining  data  for  use  in  constitutive  modeling.  This  is  mainly  due  to  the 
relative  ease  of  making  specialized  predictions  (uniaxial  strain,  plane  strain,  etc.)  from  a 
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(a)  Isotrofric  s(^. 


(b)  Transverse-isotrqnc  soil. 


Figure  7.  laical  bdiavior  of  soil  undn  hydrostatic  loading 
and  unloading  (from  Ref.  3). 
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genenlized  model  rather  than  vke  versa  (Ref.  7).  The  fdlowing  macn»c(^  continuum 
theories,  dierefoie,  are  presented  as  being  rqnesentative  of  die  basic  categories  (domination 
theory,  incremental  ]dasticity,  and  the  Endochronic  Theory  [Ref.  currendy  used  in 
modeling.  Advantages  and  disadvantages  of  each  will  be  discussed. 

1.3.1  Deformation  Theory 

Deformation  diemy  imidies  die  existence  of  a  reladvely  sinqde  relationship  between  stress 
and  strain  tensors.  The  defmmation  themy  equations  usually  i^ipear  as  follows: 


2G(o^/e^)«^ 

(2) 

3l5(o^ye^)e4j 

(3) 

where  Sy  and  ey  are  the  stress  and  strain  deviation  tensors,  reqiectivdy.  The  elastic  shear 
modulus  is  G  and  K  is  the  elastic  bulk  modulus.  Nonlinear  behavior  can  be  modeled  by  die 
introduction  of  either  stress  or  strain-dqiendent  moduli.  Onuponents  of  die  stress  and  strain 
tensors  are  0^  and  e,,,  reflectively. 

The  Duncan  and  Chang  model  (1970,  Ref.  9)  is  an  example  of  a  modd  idiidi  uses 
defmmation  dieory.  This  model  is  based  on  the  observation  that  the  stress-strain  curves  of 
soils  tested  in  a  triaxial  chamber  under  a  variety  of  cmistant  confining  stresses  resmnble  a 
hypmbda  (Fig.  8).  This  model  replicates  mily  the  axial  stress-strain  bdiavim  of  a  material 
and  assumes  a  constant  Poisson  ratio  (/t).  Improvements  to  the  model  in  1974  (Ref.  10) 
permit  the  modeling  of  nonlinear  volumetric  strains.  Additionally,  the  model  has  been 
refined  to  incmporate  behavim  observed  upon  unloading  of  die  fiedmen.  Unloading  and 
rdoading  bdiavior  is  assumed  to  be  dastic  with  the  value  of  Young’s  modulus,  dqiending 
upon  confining  pressure.  The  advantages  of  this  modd  lie  in  its  simplidty.  The  modd 
parameters  are  easily  evaluated  from  standard  laboratory  tests  and  genmally  have  a  physical 
significance  to  the  researcher.  However,  if  die  modeled  loading  path  deviates  significandy 
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from  triaxial  compression,  ttw  jnedicted  strains  may  not  be  accurate.  Since  the  model  is 
based  on  a  genmalizatior  of  Hooke’s  law  fm*  is<Mro{»c  betu  icx’,  it  cannot  modd  anisotrof^ 
responses.  Sheai  dildan.  and  posQieak  ttrain  softening  cai.not  be  modeled,  either.  Lastly, 
the  model  can  mly  jnedict  ptepalk  response. 


1.3.2  Ii 


J  Plasticity 


A  considerable  number  of  plasticity-based  models  has  been  developed  over  the  past  30  yr. 
Common  variations  include  perfect  (ideal)  plasticity,  hypodasticity,  and  viscofdastidty.  The 
most  promising  nrodds  appear  to  be  based  on  the  incremental  theory  of  plastidty  or 
viscoplastidty  in  which  the  total  strain  increment  is  expressed  as  the  sum  of  the  ela^  and 
plastic  strain  increments. 


In  most  mo(Ms,  eadi  of  these  increments  is  rate  indq)endent  (Ref.  1 1).  Regardless  of  die 
variati(m  being  considered,  the  basic  equaticm  of  all  incremental  plastidty  modds  is  as 
follows: 


dt 


(5) 


where  J,  is  the  first  invariant  of  the  stress  tensor,  K  is  the  dastic  bulk  modulus,  Sy  is  the 
stress  deviation  tmisor,  G  is  the  elastic  shear  modulus,  dX  is  coeffident  of  proportionality, 
and  f  is  die  yidd  functicm.  One  must  consider  the  ramifications  of  having  a  flow  rule  in 
Equation  S.  It  is  die  flow  rule  which  provides  for  the  calculation  of  the  plastic  strain 
increment  (de^y)  and  spedfies  whether  flow  is  associated  or  nonassociated.  The  plastic 
strain  increment  is  determined  as  follows: 


dOu 


(6) 


13 


PL-TR~92-1064 

wbefe  g  rqiresents  the  plastic  potential  function.  This  function  may  or  may  not  coincide 
with  the  yield  function  (f)  at  which  plastic  strains  develop  (Fig.  9).  In  the  case  of  associated 
flow,  the  i^astic  strain  increments  are  nOTmal  to  the  yield  surface  (i.e.,  f  g). 

Furthermore,  Drucker’s  postulate  for  stability  (1950,  Ref.  12)  aiq)lies.  However,  soil 
reqxmse  is  typically  unstable  as  the  soil  approaches  yield  conditicms  due  to  strain  softening 
bdtavior.  Thus,  the  use  of  an  flow  rule  precludes  the  accurate  modeling  of  some 

bdiavior  in  certain  soils.  Furthermore,  usage  of  the  associated  flow  rule  for 
sdls  which  exhibit  strain  softening  results  in  exc^vely  large  volumetric  strain  jnedicticms 
when  compered  to  experimental  results.  This  is  because  accurate  nuxleling  of  soils  which 
are  expected  to  uiuiergo  strain  softming  a^ypears  to  require  the  use  of  a  nonassodated  flow 
rule  (Refs.  13-15). 

Proper  modeling  of  soils  which  are  mcpected  to  strain  harden  dictates  die  inclusion  of  a 
strain  hardening  law.  The  strain  hardoiing  law  defines  modifications  to  the  yield  surface 
during  the  course  of  plastic  flow.  Two  general  types  of  hardening  laws  (isotropic  and 
kinematic)  are  suggested.  Isotrc^c  hardening  laws  dictate  that  the  yidd  surface(s) 
isotropically  expand  or  contract  about  the  hydrostatic  axis  in  principal  stress  qiace  (Fig.  10). 
This  type  of  law  was  prqiosed  for  models  developed  to  i»edict  soil  resptmse  to  monotonic 
loading.  Kinematic  hardening  laws  were  developed  to  predict  the  response  of  soils 
subjected  to  cyclic  loading.  These  hardening  laws  dictate  that  the  yield  surfu:e(s)  translate 
in  stress  space,  sometimes  changing  sh^  or  size  during  plastic  flow  (Fig.  11).  These 
stress-strain  history-induced  modifications  determine  the  constitutive  refuse  during  load 
reversals.  Modds  with  kinematic  strain  hardening  usually  agree  more  closely  with 
experimoital  results  than  modds  using  isotropic  strain  hardening  (Ref.  11). 

The  "cap  modd*  proposed  by  DiMaggio  and  Sandin’  (Ref.  16)  and  improved  by  Baladi  and 
Rdiani  (Ref.  3)  will  now  be  presented  as  an  example  of  incremental  plasticity.  This  modd 
is  based  on  a  classical  incremental  elastoplastic  strdn-hardening  theory.  It  includes  a 
stationary  ultimate  strength  (failure  surface)  which  serves  as  the  bounding  yield  surface  and 
a  strain-harctening  moving  cap  situated  between  the  ultimate  strength  surface  and  the 
hydrostatic  stress  axis  (Fig.  10).  Stress  states  on  one  instantaneous  surface,  therefore,  are 
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Figure  9.  The  plastic  potential  functi<m  (g)  a^roaches  the  yield  function  (0  as 
macroscopic  strain  increases. 
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dia.  .Jierized  by  a  common  single  value  of  a  strain  parameter  such  as  volunmtric  pla^ 
strain.  This  parameter  is  also  used  as  the  state  paramet  n*  in  the  hardening  rule.  Thus,  diis 
hardoiing  paranM^  is  the  c  racteristic  which  determiiws  die  movement  of  the  cap  surface 
accending  to  a  semimni^caliy  developed  hardening  rule  (Ref.  7).  In  an  attmnpt  to  improve 
the  model’s  ability  to  account  for  path-dq)endent  effects  and  other  features  of  anisotro|»c 
hardening  plasticity,  Reference  3  proposed  a  means  of  dtanging  the  ratio  b^een  the  major 
and  minor  ellipse  axes.  At  the  same  time,  Sandler  and  Baron  (Ref.  17)  introduced 
kinematic  hardoiing  to  enable  the  prediction  of  hysteresis  during  xtpcaied  loading  and 
unloading.  Hius,  the  plasticity  model  has  been  refined  to  include  many  realistic  aq)ects  of 
soil  behavior;  e.g.,  dilatancy,  hysteresis,  and  dtmsificatimi  untter  isotrqnc  loading. 

However,  as  previously  mentioned,  volumetric  strain  controls  the  strain-hardening  and 
volume  change  characteristics  as  of^sed  to  void  ratio  terms.  Typically,  the  void  ratio 
terms  are  the  state  variables  in  critical  state  models  developed  expressly  for  soils.  Although 
the  relati(Hiship  between  volum^c  strain  and  changes  in  the  void  ratio  is  direct,  the  use  of 
the  volumetric  strain  as  a  state  variable  instead  of  void  ratio  makes  detomination  of  the 
model  parameters  through  laboratory  testing  more  abstract. 

As  mentimied  earlior,  the  cap  model  has  limitatitms  in  the  accurate  modeling  of  cyclic 
loading  and  the  prediction  of  stress  and  strain  in  anisotropic  soils.  This  limitadtm  aq)plies  to 
othm'  models  also  (Ref.  3),  and  it  is  apparent  that  predictions  matte  using  the  cap  nKxtel  are 
most  accurate  when  based  on  isotrc^c  soils  aiul  monotonic  loading  condititms.  Several 
models  have  since  been  proposed  which  utilize  a  combination  of  isotropic  and  Idnmnatic 
hardmiing  laws. 

1.3.3  The  Endochronic  Theory 

The  material  responses  demmistrate  the  mustence  of  rate  effects  which  are  necessarily  time 
dependent.  Mottels  which  predict  the  influence  of  rate  effects  on  nrnteiial  reqxmse  usually 
invoke  an  experimmitally  determined  empirical  factor  to  adjust  the  respemse  accordingly.  A 
distinctly  differait  approach  which  considers  tinm  to  be  an  intrinsic  material  prop^  is  tte 
Endodircmic  Theory  (Valanis  and  Read,  Ref.  8).  This  theory  is  founded  on  the  coneqH  of 
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intrinsic  time.  According  to  Valanis,  time  is  a  comm<mly  used  measure  of  change.  Using 
the  example  of  a  pendulum. 


"...A  change  in  the  position  of  a  poidulum  from  one  extreme  to  the  other 
measures  half  a  poiod.  If  the  period  of  oscillation  is  one  second,  then  half  a 
second  is  the  measure  of  change  in  the  above  positicms.  In  an  ideal  pendulum, 
the  poiod  is  invariant;  i.e.,  the  temporal  distance  between  two  evoits  (in  this 
case,  the  two  extreme  positions)  is  always  constant.  If  no  one  marics  off  on 
the  real  line  points  corresponding  to  the  ends  of  half  period  (or  periods),  these 
points  are  all  equidistant.  A  continuous  time  scale  is  now  creat^  by  assuming 
that  the  density  of  events  between  any  two  points  is,  in  fact,  constant.  As  it 
happens,  the  motion  of  the  pendulum  relative  to  its  own  time  scale  satisfies 
Newton’s  second  law  of  motion.  Consider  now  another  system  which  is  not 
cognizant  of  a  Newtonian  time  scale.  A  metal  whose  stress  response  is  strain 
rate  insensitive  is  such  a  system.  Tbe  internal  changes  that  determine  the 
stress  state  begin  at  the  onset  of  deformation  and  end  when  deformatitm  stops. 
No  further  internal  changes  will  take  plac»  unless  the  state  of  deformation 
changes.  This  is  another  instance  in  which  time  is  a  measure  of  change.  In 
this  case,  however,  if  there  is  no  internal  change,  there  is  no  time  change. 
Since  the  stress  is,  inevitably,  determined  by  the  history  of  intonal  change,  it 
is  obvious  that  this  history  cannot  be  defined  with  respect  to  the  Newtonian 
time  scale  but  by  one  that  is  determined  by  the  internal  material  changes 
brought  iboMt  by  the  deformation.  Such  a  time  scale  is  obviously  intrinsic  to 
the  material  at  hand.  Thus,  the  Endochronic  theory  advocates  that  stress  is  a 
function  of  the  history  of  plastic  strain  expressed  in  terms  of  ’intrinsic  time’ 
with  this  time  being  the  length  of  the  pla;^c  strain  path  in  plastic  strain 
space." 


This  theory  is  considered  to  be  applicable  in  the  study  of  plasticity  in  granular  media 
because  sands  are  considered  to  be  strain-rate  insensitive. 


1.3.4  Local  and  Non.ocal  Theory  Formulations 

Two  general  approaches  exist  when  modeling  the  global  response  of  granular  media,  the 
local  and  the  nonlocal  theories.  A  local  theory  is  predicated  on  the  understanding  that  the 
stress  in  a  neighborhood  is  a  function  of  the  history  of  strain  in  the  neighborhood.  A 
uniform  surtiu^e  traction  will  result  in  a  linear  deformation  gradient  within  the  specimoi  and, 
thus,  a  uniform  strain  gradient.  Other  solutions  cannot  exist  without  a  loss  of  well  posedness 
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as  will  be  explained  shortly.  However,  certain  phenomena  which  cannot  be  explained  using 
a  local  theory  (e.g.,  shear  banding  and  nonlinear  deformation  gradients)  have  been  observed 
ernerimoitally.  These  pheitomena  are  termed  localizations.  It  is  well  accepted  that  these 
.omen?  can  only  be  modeled  by  the  use  of  a  nonlocal  theory  which  can  accommodate 
nonlinear  deformation  gradients.  Brief  descriptions  of  the  local  and  nonloca.  theories  are 
presented  next. 

One  must  make  an  informed  choice  when  selecting  the  partial  differential  equation  to  model 
materials  which  demonstrate  nonlinear  gradients.  Consider  the  partial  differential  equations 
which  define  the  wave  propagation  equation  and  the  diffusion  equation.  These  partial 
differential  equations  are  hyperbolic  and  parabolic,  respectively.  Equation  7  is  the  wave 
propagation  equation  (hyperbolic)  and  Equation  8  is  the  diffusion  equation  (parabolic). 

3^11  —a  _ 

--il.  +  +  sLz.  s  C*-—  (7) 

dr*  dy*  dz*  df* 

ii.il. (8) 

ax^  dy^  Si^  dt 


If  one  expresses  the  Laplacian  (V^  in  a  standard  form,  the  equations  appear  in  a  more 
recognizable  format. 


=  c* 


3^« 

dr 


(9) 


b  — 
dt 


(10) 
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The  laptecian  equals  tiie  fdlowing: 

^  ^  ^ 
dy^*  at* 

Thus,  tte  wave  equaticMi  is  characterized  by  the  wave  propagation  q)eed  (c)  ami  the 
diffusion  equation  by  the  diffiisivity  constant  (b).  Both  of  these  equaticms  are  open-«ided  in 
that  any  positive  values  of  c  and  b  are  accej^le  f(»  use  in  the  equations.  If  the  right-hand 
side  of  either  equals  a  constant,  the  equations  reduce  to  dliptical  equaticms  like: 

^  -  CONSTANT  or^T‘  CONSTANT  (^2) 

The  ell4>tical  form  is  used  to  obtain  closed-form  solutions  for  ccmditicms  of  static 
equilibrium.  It  is  the  elliptical  form  of  these  partial  diffoential  equaticms  that  is  typically 
used  to  satisfy  equilibrium  ccmditicms  in  the  formulaticm  of  nonlocal  constitutive  modds. 

The  formulaticm  is  necessarily  opea-tadoi  and  equilibrium  cannot  be  satisfied  when  one 
uses  the  higher  order  partial  diffcoential  equaticms  (either  parabolic  or  hypmbolic)  in  an 
attempt  to  mcxlel  the  aforementioned  localizations.  Attempts  to  use  the  higher  order  partial 
differaitial  equations  to  modd  the  nonlin^ties  associated  with  Icxalization  result  in  a 
formulaticm  which  is  mesh-dq)endent.  Wdl  posedness  is  lost. 

Nonlocal  theories  are  being  investigated  in  an  attempt  to  modd  localizations.  On  the 
evidence  of  some  recent  litorature  (Refs.  18  and  19),  it  af^mars  that  nonlocal  theories  have 
the  potential  of  explaining  the  cmset  of  localization.  The  results,  however,  are  not  quite 
definitive  yet.  There  is  currmtly  no  robust  theory  which  will  predict  the  onset  of  shear 
banding  without  the  inserticm  of,  for  example,  a  weakened  element  within  a  finite  demmit 
mcxiel  or  some  other  usm'-input  catalyst. 

A  qualitative  mcplanation  of  local  theory  is  provided  next  using  the  example  of  uniaxial 
loading.  Refer  to  Figure  12.  A  uniform  load  is  applied  at  the  top  of  the  specimen  and 
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Figure  12.  Linear  displacement  gradient  and  uniform  strain  prediction  of  local  theory  formulations. 
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some  vestical  di^dacement  (x)  results.  The  resulting  displacement  (U)  within  the  spedmoi 
is  thoi  a  linear  function  of  x.  They  are  directly  related  by  the  following: 

lf-p(x)  (13) 

e_  ■  ■  a(consUttit)  (14) 

*  a* 


And  thus. 


p  -  «  (15) 

Strain  (ej  at  any  point  in  the  qiecimen  is  the  partial  derivative  of  U  with  respect  to  x  and  is 
constant  throughout.  Thus,  the  local  theory  will  predict  a  linear  relationship  between 
displacement  of  the  specimen  surface  and  resulting  displacement  at  some  point  within  the 
^ledmen  and  ccmstant  strain  throughout.  Furtiwrmore,  unloading  will  result  in  a  line  which 
retraces  the  loading  path  back  towards  the  origin.  This  relationdiip  will  be  true  whether  one 
is  modeling  an  elastic,  plastic,  elastic'-visccq)lastic  material,  etc.  A  key  dement  of  this 
discussion  is  the  assumption  of  no  friction  at  the  sfpedmen  container  sidewall.  However, 
friction  is  ever  present  in  the  systems  (soil,  structure,  soil-stnicture)  whidi  these  modds  are 
developed  for,  and  to  some  extent  present  in  the  laboratory  devices  used  to  determine 
material  parameters.  Therefore,  one  cannot  neglect  its  influence.  Figure  13  shows  the 
same  uniaxial  test  chamber  of  Figure  12,  but  the  effects  of  sidewall  friction  have  been 
inferred. 

(Qualitatively,  this  influence  results  in  a  drag  force  bdng  exerted  on  the  specimen  as  it 
moves  downward  in  response  to  loading  cap  displacement.  The  magnitude  of  the  drag  force 
is  a  function  of  the  friction  coeffident  at  the  interface  as  well  as  the  lateral  force  exoted  by 
the  spedmoi  on  the  container  wall.  The  significance  of  this  influence  is  shown  in  the 
extreme  case  of  the  development  of  a  shear  band  in  granular  soils.  The  width  of  the  band  is 
typically  6-10  grain  diameters.  Within  the  band,  grains  displace,  rotate,  and  fracture, 
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Figure  13.  Nonlinear  displacement  gradient  and  strain  prediction  of  nonlocal  ttmoty  formulations. 
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resulting  in  oieigy  dissipation.  In  the  extreme,  if  one  considers  the  sidewall  interface  as  the 
middle  of  a  shear  band,  the  influence  of  the  wall  may  only  extend  into  the  ^wcimen  fw  a 
distance  of  three  to  five  grain  diameters.  As  tte  localized  disturbance  associated  with  the 
formation  of  a  shear  band  is  not  observed  at  the  sidewall  while  the  q)ecimen  is  under 
uniaxial  strain,  the  influence  of  sidewall  friction  is  thought  to  be  small.  One  must  then 
(XHisider  how  the  presence  of  the  wail  alters  the  fabric  of  the  soil  near  tl^  interface  and 
whether  this  influence  may  alter  the  perceived  micromotions.  As  discussed  in  Secticm  6.0, 
direct  comparison  of  the  average  grain  displacements  on  the  face  of  the  ^)ecimen  with  those 
of  grains  within  the  specimen  provides  an  indication  of  the  influence  sidewall  friction  has  on 
the  observed  micromotions. 

1.3.5  Summary  of  Phenomenological  Theory 

The  loading  re^nse  of  granular  material  has  been  shown  to  be  dq)endent  on  such  aspects 
as  stress  state,  loading  rate,  initial  macroscopic  conditions,  loading  path,  d^ree  of 
saturation,  and  drainage  conditions,  among  others.  Attempts  to  develop  a  model  accurately 
accounting  for  all  of  these  variables  have  frustrated  due  to  a  lack  of  understanding  of  how 
each  of  these  aspects  individually  influence  the  material  response  and  how  they  intoact  with 
one  another  as  the  original  structure  of  soil  changes  during  loading.  Therdore,  a  wide 
range  of  innovative  models  has  been  proposed  to  model  specific  aq)ects  of  the  req;Kmse  or 
predict  behavior  resulting  from  a  ^)ecified  loading  r^ime.  These  models  are  based  on  the 
principles  of  elasticity,  plasticity,  elastic  viscoplasticity,  etc.  The  variety  is  increased  by  die 
need  to  invoke  either  an  associated  flow  rule  or  a  nonassociated  flow  rule  and  the  existence 
of  localizations.  However,  all  of  these  models  can  be  coUectively  grouped  under  the 
heading  of  phenommological  models.  Such  models  attempt  to  predict  the  stress-strain 
response  of  soils  following  the  determination  of  material  parameters  which  are  subsequendy 
used  to  characterize  the  model.  These  models  frequendy  incorporate  experimentally 
determined  empirical  factors  as  a  basis  for  predicting  material  response.  At  this  writing, 
experimoits  are  underway  at  the  U.S.  Army  Corps  of  Engineers  Waterways  Experiment 
Station,  Vicksburg,  Mississippi,  to  develop  a  constitutive  law  based  on  triaxial  laboratory 
experimoits  conducted  expressly  with  the  intention  of  investigating  hardening  rules,  changes 
in  the  shape  of  the  yield  surface,  normality,  etc.  This  effort  involves  the  conduct  of 
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experimems  on  many  Ubaratoiy  ^ledmens  to  b^ter  define  a  single  yield  suriiice.  Tlie 
mults,  thou^,  will  be  ^>|dkable  only  to  amis  ^i^iich  will  be  sid>jected  to  monotonic 
.  4ding.  Thus,  die  resulting  model  may  not  accurately  predict  si^  reqxmae  during  cydu 
loading  due  to  the  dfocts  amociated  with  prestraining.  The  resulting  modd  will  be 
devdoped  baaed  on  phenmnenological  results  and  not  micrmnedianical  considerations.  The 
microstiuctural  attributes  of  the  material  being  tested  are  not  considered  in  the  formulation 
of  the  modd. 

In  summary,  some  local  theory  phenomendogical  models  have  proven  extrandy  successful 
in  predicting  soil  response  to  qiecific  loading  r^imes.  Certain  aqiects  of  soil  b^vior  sudi 
as  shear  banding  and  nonlinear  deformation  gradients  (called  localizations),  however,  cannot 
yet  be  accurately  modeled  using  such  techniques.  Nonlocal  theories  have  been  recently 
devdoped  to  accmnmodate  die  devdopment  of  localizations  and  smne  show  reasonable 
promiae.  An  alternative  iqiproadi  to  modeling  the  continuum  response  of  soils  is  die 
micnmiedianical  m  mkrostructural  qiproach.  These  models  attmnpt  to  rqilicate  die  stress- 
strain  response  of  soUs  using  dieories  founded  on  microstiuctural  aqiects  of  sml  qpedmens. 
For  granular  soils,  such  aspects  indude  grain  geometry,  tebrie,  intergranular  contact  forces, 
etc.  Basic  considerations  of  sudi  diemies  will  be  presented  in  Subsection  1.4. 

1.4  MICROSTRUCTURAL  fPARTlCULATEl  CONTINUUM  THEORY 

Microstructure-based  continuum  dieory  is  predicated  on  the  understanding  dud  die 
constitutive  response  of  a  granular  material  is  influenced  by  microstructuial  properties 
(fidnic,  interpardde  contact  laws,  granule  oonqwsition,  and  cementation,  etc.)  as  wdl  as 
mamosoopic  properties  sudi  as  void  ratio,  unit  wdght,  d^iee  of  saturation,  etc.  Sudi  an 
iqiproadi  to  modding  typically  requires  that  sini]difying  assumptions  be  made  r^arding  the 
gemnetry  of  die  granular  media.  The  media  is  usually  rqnesented  by  an  assemblage  of 
equal-sized  disks  (qdieres)  or  mixtures  of  two  or  three  different  sizes.  Furthermme,  die 
material  should  be  statistically  isotropic  with  regard  to  contad  mientadon.  These 
assumptions  are  made  to  sinqilify  the  computational  demands  of  the  reflective  modds. 
Recoit  advartees  indude  die  incoipmration  of  elliptical  (2D)  disks  or  ovoid  (3D)  shsqies 
(volumes)  in  discrete  dmnent  formulaticms.  Chang,  Misra,  and  Sundaram  (Ref.  20)  have 
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investigated  the  influence  of  miciostnictural  a;q)ects  on  the  modding  of  granular  media 
under  low  amjditude  oscillations  (<0.01%  strain).  Cundall  and  Strack  (Refs.  21  and  24), 
Pdralds  and  Dobry  (Ref.  11),  ami  Trent  (Refs.  23-25),  among  others,  have  investigated  the 
evrdutkxi  of  microstructure  in  two  and  three  dimenaons  during  arbitrary  loading  of 
computer-generated  granular  media  duough  the  use  of  discrde  dement  computm  programs. 

The  following  discussion  presorts  basic  microaxuctural  considerations  when  one 
contnnplates  modeling  the  elastic  re^ronse  of  granular  media  (wherein  no  {dastic  strain  is 
considered  to  occur)  and  modding  the  plastic  req^onse  wherein  dre  full  array  of 
microstructuial  effects  (grain  diqdaoement,  rotation,  fracture,  and  disint^ration)  are  known 
to  occur. 

1.4.1  Hastic  Response  of  Granular  Media 

A  "^ffhess  tensor*  is  commonly  used  to  rdate  stress  and  strain  in  the  generalized  Hooke’s 
law  for  elastic  deformation.  This  tensor  is  strongly  dqrendent  on  the  packing  structure  of 
crystalline  and  granular  materials  having  a  high  decree  of  symmetry.  Examples  of 
symmetric  structures  for  which  stress  and  strain  ai^iear  to  be  dqrendent  upon  the  stiffness 
tensor  are  die  simple  cubic,  fooe-centered  cubic,  and  close-packed  hexagonal  structures. 
Since  these  packings  are  not  assumed  to  rearrange  themsdves  under  dastic  deformation,  die 
material  reqxmse  is  considered  to  be  related  to  die  stiffness  tensor  of  the  symmetric 
structure.  The  same  relationship  between  stiffoess  tensor  and  material  reqxMise  is  also 
thought  to  exist  for  random  packings  of  grains  found  in  sands.  Chang  and  Misra  (Ref.  26) 
ctmclude  that  for  soils,  the  r^ularity  assodated  widi  symmetric  structure  does  not  exist  and 
"...it  is  easio*  to  derive  the  stress-strain  relaticxiship  based  on  a  microstructural  continuum 
approach.”  Therefore,  derived  stiffoess  tensors  ^^rear  suitable  (xily  for  cases  idiere  very 
small  changes  in  the  packing  structure  are  expected  during  loading.  To  accuratdy  modd 
deformations  beyond  the  elastic  limit  of  the  material,  it  is  necessary  to  understand  the 
evoluticHi  of  the  stiffness  tensor  as  loading  progresses.  Discrete  dement  techniques  are 
sqiprqpriate  alternatives  for  modeling  this  evolution.  Basic  compcMwnts  of  diserde  dement 
formulations  with  special  mnphasis  on  the  mod^g  of  frictional  materials  are  presented 
nmit. 
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1.4.2  Mastic  Response  of  Granular  Media 

Discrete  elemeiit  techniqiies  are  not  widely  acoq)ted  as  being  suitable  for  the  modeling  of 
interaction  between  discrete  bodies.  Originally  developed  fcu’  the  study  of  rock  block 
movements,  this  tedmique  has  been  refined  to  permit  the  study  of  interaction  between  lOQs 
(and  even  l(X)0s)  of  arbitrary  shapes  (vdumes)  under  a  wide  variety  of  static  and  dynamic 
loading  configurations.  The  tedinique  is  particulariy  well  suited  for  the  study  of  fiktional 
materials. 

The  principal  difference  bdween  discrete  elemoit  modding  and  classical  continuum 
modeling  is  diat  the  discrete  element  method  considers  the  specimen  to  be  made  of  10s, 
100s,  or  ICXXh  of  individual  grains  which  interact  widi  one  aiK>dier  during  loading.  This 
interaction  allows  for  highly  nonlocalized  reqxmse  in  smne  portions  of  the  ^[tecimen.  As 
mentioned  earlier,  classical  Oocal)  continuum  modds  are  typically  based  on  the  assumption 
that  uniform  macroscopic  loading/defoimation  of  a  qtedmen  results  in  a  liitear  d^ormation 
gradient  and,  thus,  a  unifimn  strain  gradient  within  the  q)ecimen.  Other  sdutions  cannot 
exist  without  loss  of  die  well  posedness  of  the  boundary  value  problem.  Thus,  discrete 
element  modds  have  been  devdoped  in  part  to  account  fim  the  possibility  that  uniform 
surfice  loadings  may  result  in  nonunifmm  interior  deformation  and  strain  gradient  fields. 
Two  such  modds,  TRUBAL  and  CONBAL,  will  be  {wesented. 

1.4.2. 1  TRUBAL.  Discrete  (di^ct)  element  programs  have  been  used  to  exercise 
constitutive  models  which  are  based  on  the  microstructural  changes  observed  in  granular 
media  under  arbitrary  loading.  The  first  widdy  accqNed  program  of  this  type  was  BALL 
(CundaU  and  Struck,  Refs.  27  and  28).  The  program  was  intended  to  account  for  the 
changes  in  fabric  which  occur  in  granular  media  under  deformation.  BALL  is  used  for 
modding  two-dimensicmal  grain  structures;  tiie  program  TRUBAL  is  used  for  modeling  in 
three  dimensions.  Many  improvements  have  beat  made  to  the  basic  program  as  well  as  the 
constitutive  model  within  it.  Such  improvements  include  the  capability  of  modding  the 
reqxxise  of  assemblages  of  different  size  grains  and  modeling  the  response  of  lightly 
cemented  grains.  Improving  the  technique  caused  modification  of  the  contact  laws  which 
dictate  the  force-displacement  relationship  between  individual  grains.  Variations  of  these 
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in^MDWBd  modds  have  been  created  to  focus  m  qwdfic  aq)ects  of  granular  media  response 
as  wdi  as  produce  ou^t  in  a  taih»ed  format.  This  subsection  will  focus  on  die 
devdopment  oi  the  original  concept  of  the  BALL  jMX^ram  and  its  inqdementation. 
Limitations  in  diis  modding  technique  will  also  be  discussed. 

The  BALL  program  uses  mi  exididt  dme-maidiing  schrnne  to  track  die  motion  of  individud 
parddes,  with  the  motion  of  each  being  governed  by  ccmservadon  of  momentum.  Figure  14 
shows  two  partictes,  x  and  y,  idiich  are  in  contact  at  point  C.  Two  unit  vectws  (q  and  tj 
are  ddined  as  pdndng  foom  the  center  of  disk  x  to  the  center  of  disk  y  and  rdated  9B' 
d^rees  dockwise,  reqiecdvdy.  The  rdadve  vdodties,  Xi,  of  die  points  cm  overisqiping 
surfinces  are  shown  to  be  functions  of  the  absolute  rotations  and  vdocities  of  each  particle. 
The  nmrmd  and  shear  cmnponents  of  these  vdodties  may  then  be  calculated  as  the 
prpiections  of  Xj  in  the  ap|gO|ffiate  directions,  e.g.,  and  reqiectivdy.  The  vdodties 
are  assumed  constant  over  a  time  stq>.  Thus,  the  incmnental  diylacements  are  obtained  by 
time  int^ration.  The  normal  and  shear  stiffnesses  (k.  and  kj  are  then  used  to  obtain 
incremental  forces.  These  are  then  added  to  die  previous  values  to  obtain  totd  nmnial  and 
shear  forces  with  the  constraint  dut  the  shear  force  is  limited  to  a  ffadion  of  the  normal 
force,  tan^„  [dus  a  constant  cdieaon.  This  sequence  of  calculations  is  performed  for  eadi 
contact  prior  to  the  qiplicdion  of  die  motkm  equations  f(v  any  partide.  Since  all  forces  are 
known  at  every  contact  point,  the  mcHnent  and  fnce  totals  on  each  particte  may  dwn  be 
used  to  calculate  translationd  and  rotational  accelerations  dirough  a  difference 
qiproximation  to  Newton’s  second  law.  Tlw  nitmients  and  fwces  are  assumed  to  be 
constant  from  time  to  so  the  vdodties,  centered  at  tinw  are  obtained  by 
int^rating  the  acoderations  over  at  and  incrementing  the  previous  values.  A  second 
integration  yidds  the  new  position  and  rmation,  centered  at  The  time  stq>  is  taken  to 
be  small  enough  diat  particles  are  unaffected  by  motions  more  than  one  diameter  away.  The 
new  total  vdodties  will  dien  be  used  to  calculate  accderations  at  time  after  the  force- 
diqdaconent  idations  have  again  been  applied  to  the  contacts.  As  inteiparticle  overlaps 
diange,  contacts  are  ddded  and  new  ones  added.  For  additional  ddail,  see  Cundall  and 
Strack  (Refs.  27  and  28). 
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1. 4.2.2  CONBAL.  Because  of  shortcomings  in  the  early  version  of  the  discrete  elemmt 
program  TRUBAL,  particularly  with  respect  to  the  interparticle  contact  laws,  P^ralds  and 
Dobry  (Ref.  11)  refined  the  code  by  introducing  Hertz-Mindlin  ctmtact  laws  (Refs.  29-31). 
This  research  was  proposed  to  develop  a  constitutive  model  for  granular  soils  expressed  in 
traditional  ctmtinuum  mechanics  terms  but  grounded  on  the  particulate  micromechanical 
behavior  of  the  material.  The  primarily  analytical  effort  involved  refinement  of  the  laws 
modeling  the  response  between  two  elastic  spheres,  using  the  new  laws  to  study  the  response 
of  regular  and  random  arrays  of  spheres.  The  refined  contact  laws  were  incorporated  in  a 
discrete  element  program.  Analytical  predictions  using  the  program  were  compared  with 
test  results  on  glass  beads  subjected  to  axial-torsimial  tests  in  a  hollow  cylindrical  loading 
device.  The  experiments  and  numerical  simulations  confirmed  that  the  yield  surface(s)  of 
the  granular  material  compress(es)  along  the  loading  axis  as  loading  progresses  but  not  in 
other  directions. 

As  mentioned  previously,  this  effort  required  generalization  of  the  contact  laws  for 
interaction  between  two  rough  spheres.  Although  Hertz,  Mindlin,  and  Deresiewicz  (Refs. 
29-33)  had  solved  this  problem  for  a  limited  number  of  cases,  incorporation  of  the  law  into 
a  discrete  element  program  required  laws  to  be  developed  for  arbitrary  loading  and  the 
utilization  of  an  incremental  elastic-plastic  model. 

The  question  of  contact  law  between  two  particles  is  most  critical  in  modeling  the 
deformation  of  noncemented  granular  materials.  This  is  because  the  macroscopic  re^Kmse 
is  considered  to  be  a  function  of  the  elasticity  of  the  particles  and  the  geometry  of  the 
contacts  as  well  as  the  coefficient  of  friction  at  the  contact.  Particle  elasticity  and  contact 
geometry  appear  to  be  the  most  important  at  low  levels  of  strain  with  particle  sh2q)e,  contact 
geometry,  and  friction  becoming  the  controlling  factors  as  interparticle  shearing  increases. 
Reference  29  demonstrates  that  the  normal  force  deformation  behavior  at  a  contact  is 
nonlinear  elastic  if  fracture  is  not  permitted.  Various  solutions  have  since  been  proposed  for 
contact  between  unequal  diameter  spheres  with  arbitrary  combinations  of  normal  and  shear 
forces  being  applied.  Typically,  these  solutions  dictate  an  infinite  shear  stress  at  the  edge  of 
the  contact  area  and  permanent  slip  is  induced.  This  aspect  of  the  model  introduces  the 
need  to  consider  energy  dissipation  and  represents  the  key  advance  in  Seridi  and  Dd)ry’s 
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sdutkm  (Ref.  34)  vice  the  early  modds  of  Mindlin  (Refs.  30  and  31).  Figure  IS  shows  the 
nonlinear  fbfce-defonnation  curve  associated  with  the  applicatimi  of  a  monotonically 
increasing  tangential  force  to  two  q)heres  in  contact  under  a  constant  normal  force. 
Reference  34  luis  solved  this  problem  with  the  nmstitutive  law  founded  on  the  incremental 
themy  of  plasticity.  The  model  assumes  kinematic  hardoiing  and  relies  on  the  existence  of 
an  infinite  number  of  yield  surfaces.  These  are  features  which  have  been  shown  to  jxovide 
the  most  accurate  predictions  in  eariier  phenomoiological  models.  The  ccmstitutive  rdatitm 
(implemented  in  the  program  CONTACT)  has  been  shown  to  provide  accurate  rq>roductions 
of  the  analytical  solutions  obtained  by  Mindlin  and  D^edewicz  (Ref.  31).  The  constitutive 
law  has  since  been  incorporated  in  the  discr^  element  program  TRUBAL  resulting  in  die 
program  CONBAL. 

Altiiough  the  aforementioned  modifications  in  the  discrete  elemoit  program  BALL  have 
resulted  in  improved  predictive  oqnbilities,  it  appears  that  a  wider  variety  of  nucro- 
structural  effects  must  be  incorporated  before  one  can  predict  respcmses  in  a  complex 
loading  oivironmoit.  As  a  minimum,  localized  displacement  of  particles,  the  development 
of  clusters,  individual  particle  fracture  and  disintegration,  and  particle  rotation  must  be 
explored  and  incorporated.  If  moist  or  saturated  soils  are  to  be  modeled,  then  one  must  also 
OHisider  effective  stress,  hydrostatic  pressure,  dynamically  included  pore  pressure,  and  fluid 
flow  within  the  voids.  Many  attraipts  have  been  made  to  quantify  changes  in  these 
nucrostructural  parameters  on  specimens  whidi  have  beoi  subjected  to  some  pre^pecified 
macroscopic  strain  or  stress  level.  Some  of  these  efforts  have  evoi  provided  such 
information  as  the  soil  (soil  simulant)  was  deforming.  However,  limitations  in  mcperimratal 
techiuques  and  computational  abilities  have  precluded  the  development  of  a  unified  theory 
which  links  microstructural  effects  to  the  continuum  re^nse  of  frictional  matoials. 

1.4.3  Summary  of  Microstructural  Theory 

Sevoal  computer-based  techniques  have  been  used  to  aid  in  the  investigatimi  of 
microstructural  effects  and  their  influoice  on  the  continuum  resprmse  of  granular  media. 

The  principal  techiuque  involves  using  the  discrete  element  codes  BALL  and  TRUBAL  (or 
their  derivatives)  to  conduct  parametric  studies  on  macroscopically  similar  qredmens. 
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Elastic  sirfieres  under  normal  and  tangential  loads. 


Figure  IS.  Nonlinear  force-defonnation  curve  for  elastic  spheres 
under  normal  (N)  and  tangential  (T)  loads  (constant  N, 
T  increasing)  (from  Ref.  34). 
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Similarly,  numerical  studies  are  conducted  on  "specimens*  which  have  been  changed  subtly 
or  siri)jected  to  various  stress  paths  prior  to  testing.  The  purpose  of  such  testing  is  to 
identify  and  quantify  the  influence  of  some  specified  microstructural  effect.  A  perceived 
limitation  is  the  inability  to  model  other  than  disks  (in  2D)  and  ^heres  (in  3D).  Thus,  the 
infltmice  of  particle  stuqx  could  not  be  investigated  using  the  early  codes.  Recent  advances 
in  the  discrete  dement  technique  involve  the  incorporation  of  particle  bmiding  and  the 
ability  to  modd  elliptically  shsqped  particles.  These  codes  provide  the  stress*strain  roqnmse 
of  the  assemblage  as  well  as  the  response  to  elastic  wave  passage;  however,  only  a  limited 
aqpability  exists  to  highlight  and  visualize  specific  microstructural  evoits  which  may  control 
the  macroscopic  response.  As  the  codes  are  refined  to  include  irregularly  sh£q)ed  particles 
and  particle  fiacture,  the  ability  to  isolate  and  study  ^)ecific  microstructural  occurrences 
may  provide  insight  n^arding  the  formation  of  localizadtuis  which  control  the  overall 
response. 

1.5  £12B£QS£ 

The  purpose  of  this  research  was  to  investigate  the  evolution  of  microstructural  effects  in 
granular  material  under  uniaxial  strain  of  up  to  10.0  percent.  The  relative  influence  of 
several  microstructural  effects  (such  as  grain  di^lacement,  rotation,  fiacture,  and 
disintegration)  on  strength  mobilization  was  studied.  These  microstructural  effects  were 
investigated  to  quantify  the  influence  of  each  on  the  constitutive  response  of  graiiular 
material  and  to  devdop  methods  of  incorporating  such  effects  into  a  microstructure-based 
constitutive  model.  Special  attention  was  paid  to  those  microstructural  effects  which  easily 
lend  themselves  to  incorporation  in  a  discrete  element  computer  code. 

1-6  OBJECTIVES  AND  SCOPE 

The  principal  objective  was  to  investigate  the  evolution  of  fabric,  the  preferred  orientation 
of  individual  grains  as  well  as  their  packing  arrangement,  in  granular  media  under  uniaxial 
strain.  T?  purpose  was  also  to  investigate  the  evolution  of  damage  in  granular  media  under 
uniaxial  strain  to  better  understand  how  these  microstructural  effects  contribute  to  the 
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ccmstitutive  re^xnise  of  the  material.  A  related  goal  was  to  quantify  microstructuial  effects 
which  occur  on  the  surface  of  a  granular  specimen  as  weU  as  intonal  effects.  Such 
information  is  required  to  discmn  how  the  exigence  of  a  q)ecimen  txnmdary  might  influence 
the  manifestation  of  microstructuial  effects.  Furthermore,  this  effort  was  intended  to 
identify  the  :nicrostructural  effects  which  most  directly  influence  the  macroscq>ic  reqxmse 
of  granular  media.  Lastly,  a  goal  was  to  idoitify  methods  of  refining  studies  of  the 
evolution  of  fltbric  in  granular  media  under  loading. 

Cuboid  specimens  of  poorly  graded,  well-rounded  silica  grains  (nominal  diameters  of  O.SO 
and  0.7S  mm)  were  subjected  to  uniaxial  strain  in  a  rigid-walled  confinement  vessel. 
Micromotions  and  the  evolution  of  damage  were  directly  obsoved  through  sapphire 
windows  in  the  sides  of  the  confinement  vessel  in  a  series  of  displacement-controlled  tests. 
Companion  displacement-controlled  tests  were  conducted  using  x  rays  to  study  the 
micromotion  of  embedded  lead  particles.  Macroscopic  vertical  strain  was  limited  to  10 
percent.  Lateral  strain  was  not  permitted.  Selected  specimens  were  qmxy-impr^nated  and 
subsequently  sectioned  to  investigate  internal  microstructuial  changes.  The  microstiuctural 
effects  observed  during  loading  and  following  epoxy  impiegnaticm  were  related  to  the 
constitutive  response  of  the  material  and  limitaticms  in  the  current  experimental  techniques 
(principally  edge  effects)  are  discussed.  An  int^iated  mcperimental/computational 
nucrostructure  study  technique  is  proposed. 
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2.0  REVIEW  OF  PREVIOUS  RESEARCH 

Sevenl  distinct  approaches  have  beoi  taken  r^arding  the  s^y  of  the  nucrostructuial 
reqmnse  of  loaded  soils.  The  most  direct,  but  potentially  tedious  an>roach  involves 
quantifying  the  coordination  number  (numbo'  of  contacts  per  particle),  (xmtact  orientation, 
and  changes  in  these  attributes  following  load  aj^lication.  This  approach  generally  involves 
immersing  a  ^)ecimen  in  paint  or  ink  and  thoi  draining  the  pore  fluid  to  create  rings  at 
particle  contacts,  the  number  and  orientation  of  which  can  subsequently  be  quantified. 
Another  common  method  of  microstructural  investigatimi  involves  qmxy  impr^nation  of 
soil  specimois  to  permit  sectimiing  and  direct  observation  of  specimen  internal  attributes. 

Yet  another  approach  involves  the  use  of  embedded  lead  particles  and  x  rays  to  directly 
view  or  infer  the  displacement  (strain)  which  occurs  in  a  specimm  during  loading.  La^y,  a 
common  technique  for  studying  the  loading  response  of  soils  involves  the  use  of  assemblies 
of  photoelastic  disks  subjected  to  various  loading  regimes  while  being  viewed  in  polarized 
light.  This  technique  permits  the  direct  determination  of  the  stress  distribution  among  disk 
assemblies  under  either  static  or  dynamic  loading.  The  advantages  and  disadvantages  of 
these  approaches  will  be  presented. 

The  technique  of  using  paint/ink  infiltration  te^niques  to  identify  particle  ccmtacts  and  infer 
assemblage  geometry  oijoyed  brief  pqnilarity  in  the  eariy  stages  of  microstructural  studies. 
This  technique  required  the  creation  of  an  assemblage,  of  either  ^herical  or  slightly 
irregularly  shaped  particles,  which  was  subsequoitly  immersed  in  dilute  paint  or  ink. 
Drainage  of  the  pore  fluid  resulted  in  small  rings  remaining  at  the  points  of  contact  between 
adjacent  particles.  "Dissection"  of  the  q)ecimcn  permitted  the  researcher  to  determine  the 
number  of  ccmtacts  per  particle  and  the  goioal  distribution  of  contacts  on  the  particle 
surface.  The  static  nature  of  the  data  collection  technique  limited  this  type  of  analysis  to 
studies  of  the  packing  geometry  and  provided  no  information  regarding  the  evolution  of  this 
geometry,  or  "fabric,"  as  the  specimen  was  subjected  to  stress  or  ^rain.  Furthermore,  the 
contact  geometry  associated  with  each  individual  particle  served  as  the  basis  for  the 
development  of  a  geometry  which  was  then  attributed  to  an  "average"  particle.  Subtle 
changes  in  particle  orioitation  or  contact  geometry  which  might  influence  the  constitutive 
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reqKXise  of  the  assemblage  were  being  lost  during  data  reduction.  This  technique  has  also 
been  iqipUed  to  soil  specimens  (irr^ulariy  shaped  particles  with  a  nominal  diam^er  of 
~  1.0  cm);  however,  the  sheer  tedium  of  this  method  and  the  limited  data  available  make 
this  an  unattractive  technique. 

Matsuma  Oda  (Ref.  35)  is  omsidered  to  be  the  pitmeer  in  the  study  of  soil  fabric  using  the 
technique  of  epoxy  impregnation  and  sectioning.  These  studies  hdped  define  the  concept  of 
fri>ric  to  be  discussed  in  Section  3.0.  Oda’s  technique  involved  the  epoxy  impregnation  of 
soil  specimens  vriiich  served  to  "freeze"  the  individual  particle  (moitations.  Subsequent 
sectioning  of  the  ocimens  (through  the  use  of  a  diamtmd  saw)  permitted  the  researdier  to 
view  particle  orientations  and  infer  contacts.  This  technique  was  used  to  study  a  vari^  of 
granular  soils  with  literally  thousands  of  individual  particle  orientations  and  contacts  being 
observed/inferred.  Although  innovative,  several  shortcomings  in  this  technique  became 
apparent.  Polishing  the  cut  surface  of  the  q)ecimens  permitted  limited  observation  of  the 
contacts  which  existed  slightly  bdow  the  exposed  surface.  This  was  useful  as  very  few 
contacts  fell  exactly  on  the  plane  of  the  exposed  section.  Out-of-plane  contact  orientations 
had  to  be  inferred  through  manual  observation.  This  aq)ect  of  the  technique  introduced 
subjectivity.  Furthermore,  limitations  in  visualizing  the  three-dimensional  fabric  or 
incorporating  such  information  in  a  computer  code  limited  presentatitxi  of  results  to  a  two- 
dimensional  format.  Figure  16  riiows  a  rose  diagram  typically  used  in  the  iqnesentation  of 
contact  orientations.  Such  diagrams  have  become  useful  in  the  visualization  of  the  force 
exerted  on  an  "average"  particle  at  a  specified  orientation.  Thus,  Oda’s  work  highlighted 
the  need  to  explore  and  visualize  the  three-dimensional  nature  of  the  particle  contact 
orientations.  However,  it  was  limited  by  the  need  to  manually  interpret  particle  and  contact 
geometry  and  an  inability  to  visualize  fabric  in  its  true  three-dimoisicmal  form. 

Regarding  the  use  of  x-ray  techniques.  Reference  4  defined  the  strain  field  which  occurs  in 
cylindrical  specimens  of  dry,  doise  sand  subjected  to  triaxial  compression.  The  deformatitm 
field  in  laboratory-prq)ared  specimens  was  quantified  by  placemoit  of  lead  beads  within  the 
qpecimoi  in  a  vertical  plane  which  was  parallel  to  and  passing  through  the  specimen 
longitudinal  axis.  X  rays  taken  during  loading  were  used  to  quantity  the  deformation  and 
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Figure  16.  Typical  rose  diagram  for  granular  material.  The  V-secticm 

shows  that  the  long  axes  of  particles  have  assumed  a  primarily 
horizc  d  orientation  (from  Ref.  35). 
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strain  fields.  Donan  (R^.  4)  determined  that  when  deformations  are  nearly  homogeneous, 
doise,  granular  sc^s  exhibit  little  or  no  strain  softening.  This  d^ermination  was  made  by 
comparing  the  re^xxise  of  tiiaxial  qjedmens  with  some  having  either  lubricated  or 
nonlulmcated  end  platens  (Fig.  17).  Triaxial  compresaon  tests  were  poforroed  on  flat 
q)ecimens  of  dense,  dry  sand  and  concluded  that  only  an  insignificant  amount  of  strain 
softening  occurred.  Vardoulalds  (Ref.  18)  used  x-ray  radiograirfis  on  sand  q)ecimens  to 
study  deformation  pattms  in  sand  under  biaxial  compression  loading.  Macroscopically,  die 
specimm  deforms  uniformly  until  failure  is  reached;  however,  there  is  inhomogeneous 
internal  deformation  and  thin  shear  zones  develc^  and  dilate  during  loading.  The  density  of 
these  zrnies  is  qualitatively  different  than  that  of  the  remaining  sand  specimen.  CHien, 

Wang,  Schreyer,  and  Rutland  (Ref.  36)  have  used  real-time  x-ray  radiogn^hy  and  image 
analysis  to  investigate  the  growth  of  small  zones  of  localized,  rdadvely  large  deformations 
within  uniaxially  loaded  concrete  mortar  specimens.  C;onsidend>le  localized  damage  occurs 
within  qiecimens  (as  reflected  by  macroscopic  softening)  when  as  litde  as  40  percent  of  the 
peak  load  has  been  applied.  This  result  rqiplies  when  the  ends  of  the  iqiecimen  are 
lubricated  at  the  platen.  Chen  et  al.  conclude  that  the  stress-strain  relation  cannot  be 
directly  derived  fiom  the  force-di:q>lacement  curve  due  to  the  nonhomogeneous  material 
response.  This  type  of  research  has  provided  a  qualitative  glimpse  into  the  distributitm  of 
displacement  and  strain  within  laboratory-prqjared  ^)ecimens;  howevo’,  the  microstructural 
events  which  ctmtribute  to  the  development  of  the  macroscopic  response  have  ccmtinued  to 
elude  researchers. 

Photoelasticity  has  been  used  quite  effectively  in  the  study  of  particulate  matorial  response  to 
both  static  and  dynamic  load  plication.  R^arding  static  loading,  CXla,  Kishino,  and 
Nemat-Nasser,  among  many  others,  have  explored  the  response  of  granular  matoials  using 
two-dimensional  assemblages  of  photoelastic  disks.  The  disks  resided  in  a  loading  frame 
which  subjected  the  assemblage  to  a  variety  of  loading  paths  while  pomitting  the  resulting 
isochromatic  lines  to  be  studied  and  jdiotogr^hed.  These  studies  have  typically  involved 
quasi-static  loading  rates  and  circular  "particles"  of  uniform  aze.  Round  particles  of 
different  diameters,  however,  have  been  mixed  in  the  same  q)ecimen  and  dlipsoidal 
particles  have  also  been  studied.  Researdi  on  such  particle  assemblages  is  conducted  to 
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TRIAXIAL  COHPrISSlON  OF  OCKSE  SANO 


if 


•a 


Stress  ratio  versus  axial  strain  for  dense  sand  —  triaxial 
conpresslon  {after  Oenan,  1975) 


Figure  17.  Internal  defcmnation  fidd  for  dense  sand:  (a)  Lubricated  end 
I^atens,  and  (b)  Nonlubricated  end  platens  (from  Ref.  38). 
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betto'  understand  interparticle  force  and  mobilized  angle  of  friction.  Additionally,  the 
technique  pomits  continuous  observation  of  the  assemblage  leqxmse  at  the  "particulate” 
level.  The  technique  of  testing  round  photoelastic  particles  lends  itself  most  directly  to 
comparison  of  mcperimental  renilts  widi  predictions  made  using  disoete  element  programs 
such  as  TRUBAL  and  CONBAL.  But  direct  comparison  of  the  results  of  a  iriKHodastic  disk 
experiment  with  laboratory  soil  qsedmens  highlights  the  rdatively  stiff  response  of  the 
iriiotoelastic  disk  assemblage.  The  two-dimensional  rqnesentation  of  the  sml  particks  does 
not  caqtture  die  realities  of  the  three-dimensional  soil  qieciroen.  Additionally,  die 
experiments  conducted  using  photoelasdc  disks  inovide  little  or  no  information  regarding  die 
failure  and  postfailure  response  of  the  assembU^e  as  the  disks  are  typically  not  stressed  to 
the  pdnt  of  fracture.  This  technique  is  limited  by  the  need  to  conduct  tests  using  only  two- 
dimensional  rqnesentations  of  simple  geometric  shapes  and  the  inability  of  testing  beyond 
"fulure"  of  the  assemblage. 

Photoelasticity  has  been  used  for  the  study  of  wave  propagation  in  granular  media  by  Shulda 
and  Sadd  (Ref.  39).  Cinema  {diotography  using  stroboscofnc  light  permits  the  recording  of 
shock  wave  passage  through  an  assemblage  of  two-dimensional  photoelastic  disks.  As  diis 
effort  is  concerned  primarily  with  quasi-static  loading  of  particles,  microstructure  study 
using  dynamic  jdiotodastic  techniques  will  not  be  discussed. 
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3.0  Tlk.  C(»4CEPT  OF  FABk^.  IN  GRANULAR  MEDIA 

Thete  are  two  basic  ctmqxxients  c  abric,  particle  orientatkMi  and  packing  (Fig.  18),  in 
non^iberical  particles.  These  two  concqks  are  used  to  quantify  the  spatial  arrangonent  of 
partides  as  well  as  the  resulting  voids.  Within  these  two  categories,  there  are  five 
commonly  acoqMed  elmnents  of  fabric.  Particle  (mentation  can  be  quantified  wnng  the 
dementt  of  vecUM’  mean  direction  and  vectew  m?.  nitude.  Packing  can  be  quantified  using 
die  probadnlity  density  function  of  normals  at  contacts,  the  mean  of  the  coordination 
number,  and  the  standard  deviation  of  le  co(»dination  number  (Ref.  35).  These  terms  will 
be  e'  olained  sh(Mtly.  These  elmnents  are  useful  in  explimng  scune  properties  of  granular 
materials  such  as  anisotropy  and  dilatancy  during  shear.  Furdiennore,  diey  have  been  used 
as  the  basis  fm  theewies  r^arding  the  constitutive  reqxmse  of  granular  media. 


Figure  18.  ConeqH  of  fabric  in  granular  material  (from  Ref.  35). 

Particle  orientatitm  fabric  is  detomined  from  knowledge  of  the  vector  mean  directiem  and 
magnitu^  of  typical  dimoisicnial  axes  of  a  suitable  number  of  particles  within  a  given 
sample.  The  exact  determinatiem  of  the  true  long  axis  of  a  particle  is  quite  difficult  in 
practice.  Therefore,  an  alternate  procedure  has  been  developed  which  utilizes  the 
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information  available  ftom  the  thin  sectioning  of  ^xy-impr^nated  granular  iqwdmens. 
Particle  cross  sections  are  exposed  using  this  technique  and  die  long  axis  of  eadi  of  the 
exposed  cross  sections  is  designated  the  ayiparent  long  axis  of  die  particle.  The  apparent 
long  axes  are  then  summed  with  reflect  to  length  and  orientation  and  a  wdghted  avoage  of 
their  distribution  is  used  to  produce  a  rose  diagram.  The  principal  attributes  of  this  diagram 
are  the  vector  mean  directicHi  (6)  which  dmiotes  the  fneferred  directitMi  of  the  apparait  long 
axes  of  each  of  the  particles  ccmsiitered  and  the  vector  magnitude  (VM)  which  is  a  measure 
of  parallel  alignment.  Thda  and  VM  are  defined  as: 
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where  n  is  the  number  of  measurements  and  is  the  inclination  of  the  apparent  long  axis 
with  respect  to  a  specified  axis.  The  vector  ma^tude  is  expmsed  as  a  percentage  with 
0  percent  being  a  totally  random  orientation  and  100  percent  being  a  totally  parallel 
alignment.  Determination  of  these  parameters  on  three  mutually  orthogonal  planes  can  be 
used  to  define  orientation  fabric  in  three  dimoisions.  An  obvious  limitation  of  this 
technique  is  that  the  orientation  and  magnitude  of  the  apparent  long  axes  may  be  influenced 
by  the  orientation  of  the  chosen  cutting  plane(s).  Additionally,  the  tedious  nature  of 
obtaining  such  information  precludes  the  processing  of  a  large  number  of  samples. 

Packing  fabric  is  determined  from  knowledge  of  the  particle’s  positions  with  req)ect  to  one 
another  as  well  as  the  orientation  and  distribution  of  their  contacts.  Figure  19  shows  (me 
method  of  specifying  the  orientation  of  the  two  contact  normals  which  result  whoi  two 
particles  touch.  Whmi  a  substantial  number  of  particles  are  considmed,  the  ccmrdination 
number  (number  of  ccmtacts  per  particle)  becomes  an  important  paramet'^r  in  the 
determination  of  packing  fabric.  The  importance  of  the  contact  normal  orientations  and  the 
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mean  and  standard  deviation  of  dm  cocndinadon  number  are  considered  in  the  works  of 
0(fa^  Nenud-Nasser,  and  Chang  (Refs.  10,  20,  26,  35,  and  38)  anuxig  odiers,  and  will  not 
be  inesented  here.  Only  a  limited  c^nbility  exists  to  relate  static  load  bearing  capacity  to 
various  fitbrk  attributes  and  configurations.  The  microstnictural  drf(»mations  and 
micromotitms  of  individual  particles  in  granular  media  under  arbitrary  macroscopic  strain 
have  yet  to  be  incorporated  in  a  microstnicture-based  mot^.  Thus,  there  is  no  data  rdatmg 
the  evolution  of  microstnictural  effects  to  the  nuderial  constitutive  response.  This  research 
has  been  conducted  to  provide  such  information  by  the  direct  obsmvation  of  granular 
material  under  uniaxial  strain  and  also  to  provide  insight  into  dm  evcriutum  of  fabric 
the  same. 


I )  Week  ctncotrotMa  SirMtf  (•nctairatlon 

el  N|  to  <r,  el  N|  te  ffi 


Figure  19.  Schematic  representations  of  two^menaonal  granular  assemblies. 
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4.0  EQUIPMENT  DESIGN  AND  OPERATION 

This  section  (tescribes  the  equipmoit  and  some  techniques  used  to  observe  the 
microstructural  leqxNise  of  poorly  graded  sand  umto  uniaxial  macioscoinc  strain.  M^hods 
used  to  quantify  the  microstructural  effects  visible  on  the  specimen  surface  as  well  as  within 
the  q;)eciroen  are  presrated. 

4.1  SPECIMEN  CONTAINER 

A  specimen  container  (confinement  vessel)  was  designed  to  prevent  lateral  expansion  of  the 
specimen  and  permit  direct  viewing  (through  ssqqphire  windows  fixed  in  the  steel  wall)  of  a 
soil  mass*  re^nse  to  uniaxial  strain.  The  crmfinement  vessel  consisted  of  a  base  plate, 
four  walls,  and  a  loading  esq*.  The  sidewalls  extended  up  beyond  the  specimen  and  served 
as  a  guide  for  the  cap  which  was  used  to  displace  the  top  of  the  q)ecimen.  The  base  plate, 
walls,  and  loading  cap  were  made  of  4130  heat-treated  steel  to  resist  abrasion  by  the  soil 
grains.  Initial  testing  was  conducted  in  a  container  which  permitted  the  simultaneous 
recording  of  specimen  response  through  three  orthogonally  placed  sapphire  windows  (Fig. 
20).  Tests  revealed  that  no  sqipreciable  activity  occurred  at  the  base  of  any  spedmeas  up 
through  the  application  of  10  percent  macroscopic  strain.  Subsequently,  testing  was 
conducted  in  spedmoi  containers  which  permitted  simultaneous  viewing  of  the  qiedmen 
response  through  windows  in  three  of  the  sidewalls.  The  containers  were  made  to  house  a 
cuboid  ^)ecimen  with  each  side  being  50  mm  in  loigth.  The  containers  woe  designed  to 
permit  viewing  at  0.375,  0.50,  0.56,  and  0.75  specimen  height  (h)  (Table  1). 

Each  viewing  port  was  25  mm  in  diameter  and  fitted  with  a  37.5-mm-diam.  0-deg  saiq)hire 
window.  The  windows  provided  the  means  of  viewing  microstructural  effects  (individual 
grain  displacement,  rotation,  fracture,  and  disintegration)  during  loading.  Zero-d^ree 
sapphire  was  required  for  the  windows  because  this  material  did  not  introduce  any  visi^ 
distortion  as  lateral  loading  increases  on  the  window  during  uniaxial  compression  of  the 
specimen.  Fiducial  marks  were  laser-etched  on  one  side  of  each  of  the  sapphire  windows  to 
serve  as  references  agsdnst  which  grain  motion  could  be  measured.  As  the  marks  create  a 
weakness  in  the  window,  they  had  to  be  on  the  inside  (compression  side)  of  the 
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Figure  20.  Confmement  vessd  details. 
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windows  to  reduce  the  potential  for  window  ftacture.  This  was  beneficial  as  it  also  redu(^ 
parallax  error  during  micromotion  recording.  The  fiducial  mark  geometry  is  shown  in 
Figure  20.  The  center  cross  hair  in  each  window  was  located  at  the  window  crater  height 
as  listed  in  Table  1. 


Table  1.  Height  of  viewing  port  craters  for  the  specimen  containers 
window  crater  height  (mm). 


h 

0.375 

0.50 

0.56 

0.75 

50 

19.0* 

25.0 

28.0 

37.5 

NOTE:  This  window  was  set  at  19.0  mm  because  it  was  the  lowest  the  window 
could  be  set  vnthin  a  confinemrat  vessel  wall.  The  window  crater 
height  correqmnds  to  0.37S  h. 


4.2  SPECIMEN  LOADING  FRAME  AND  MACROSCOPIC  PROPERTY  RECORDINC 


All  tests  were  conducted  with  the  specimen  being  subjected  to  a  constant  macroscopic 
uniaxial  strain  rate  of  O.S  percent  per  minute  whether  loading  or  unloading.  Di^lacemrat 
control  of  the  loading  platen  was  the  chosen  method  of  ai^lying  strain  to  the  q)ecimra. 
Thus,  for  the  SO-mm-tall  specimen,  the  loading  cap  was  displaced  at  0.25  mm/min 
(0.0417  mm/s). 


Two  differrat  loading  frames  were  used  for  testing.  The  capacity  of  each  was  1(X)  kip. 
The  maximum  force  required  of  the  loading  frames  was  ~  35  kip  (for  specimens  subjected 
to  ~  12.0  percent  macroscopic  strain).  The  principal  loading  frame  was  located  at  Los 
Alamos  National  Laboratories,  Los  Alamos,  New  Mexico,  and  was  used  during  cinema 
recording  of  miraomotions.  The  second  frame  was  located  at  the  University  of  New 
Mexico,  Albuquerque,  New  Mexico,  and  was  used  when  x  rays  were  utilized  to  determine 
the  internal  micromotions  during  loading.  Standard  XY-plotters  were  used  to  display 
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measuied  load  versus  cs^  di^lacemoit  during  testing.  Addidtmally,  the  load  versus 
di^laconent  refuse  was  digitally  recorded  as  tests  progressed. 

4.3  MICROSTRUCTURAL  RESPONSE  RECORDING 

The  observation  and  recording  of  several  microstructural  effects  (individual  grain 
displacement,  rotation,  fracture,  and  crushing)  were  considered  critically  important  during 
testing.  The  distinction  is  now  made  between  the  microstructural  effects  which  are  visible 
and  can  be  recorded  during  testing  (called  surfidal  microstructural  effects)  and  the  effects 
which  can  only  be  observed  posttest  through  the  use  of  epoxy  impr^nadon  and  secdoning 
(called  internal  microstructural  effects).  Although  these  effects  are  essendally  similar  on  the 
surface  as  well  as  inside  the  specimen,  the  techniques  needed  to  observe  each  are  different 
and,  thus,  the  disdncdon  is  made.  Grain  displacemmt  and  rotadon  on  the  specimoi  surface 
were  observed  through  the  mutually  exclusive  use  of  cinema  or  video  recording  of 
micromodons  during  testing.  Internal  micromodon  was  observed  through  the  use  of  x  rays 
as  presented  in  Subsecdon  4.3.2.  Evidence  of  fracture  and  crushing  (considered  as  artifacts 
of  the  loading  process)  was  observed  after  testing  through  the  use  of  epoxy  impr^nadon 
and  secdoning  techniques,  presented  in  Subsecdon  4.3.3.  The  following  subsecdons  present 
the  two  techniques  explored  for  the  recording  of  the  surdciai  microstructural  response 
during  loading. 

4.3.1  Surficial  Microstructural  Response 

Surficial  micromodons  were  recorded  primarily  through  the  use  of  either  16-  or  3S-mm 
cinema  or  video  techniques.  The  data  regarding  grain  displacements  were  obtained  solely 
through  analysis  of  cinema  images.  Unsuccessful  attempts  were  made  to  quandfy  grain 
rotadons  from  these  images.  Shortcomings  in  this  technique  highlighted  the  need  to  develop 
a  more  refined  system.  Subsequendy,  video  recording  techniques  were  explored.  No 
quandtadve  data  were  obtained  using  the  video  techniques;  advantages  of  and  recent 
advances  in  both  techniques  will  be  presented. 
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4.3. 1.1  Data  Capture  Through  Cinema  Photography.  Cinema  images  of  the  grain 
micromotions  were  obtained  using  either  a  16-mm  system  with  (me  image  being 
recorded  per  frame  or  a  3S-mm  system  with  four  images  per  frame  recorded  simultaneously. 
The  16-mm  system  relied  on  erne  l«is  per  camera  with  one  set  of  lens  and  camera  being 
required  for  each  sapphire  window  being  viewed.  The  3S-mm  system  cemasted  of  a 
I^oto-optic  cable  which  permitted  the  simultaneous  recording  of  up  to  four  images  pa 
frame.  The  advantage  of  the  photo-optic  cable  system  (called  the  hydra-head)  was  that  the 
microstructural  response  could  be  recorded  from  several  positions  on  the  same  specimen 
with  no  ambiguity  regarding  the  application  of  macrostrain  at  the  time  the  image  was 
obtained.  Furthermore,  only  die  hydra-head  system  was  csynble  of  recording  images  from 
the  bottom  of  any  given  specimen.  One  "channel”  of  the  photo-optic  system  was  usually 
used  to  record  the  force  versus  deformation  history  from  an  XY-plotter.  Thus, 
microstructural  response  could  be  more  easily  correlated  with  macrostrain-state.  Both 
systems  utilized  the  same  SO-mm  fixed-focal  length  lenses  with  the  appropriate  extension 
tubes  for  the  desired  degree  of  image  magnification.  Although  the  hydra-head  offered  the 
advantage  of  flexibility  with  regard  to  the  choice  of  windows  being  viewed,  as  well  as  the 
ability  to  record  the  response  at  several  locations  at  once,  the  images  provided  by  the  16- 
mm  system  were  far  superior.  As  a  result,  the  hydra-head  system  was  not  used  following 
the  determination  that  there  was  no  appreciable  microstructural  activity  to  be  studied  through 
the  sapphire  window  at  the  bottom  of  the  specimen.  The  photo-optic  cable  had  an 
insufficient  number  of  strands  to  transmit  a  picture  of  quality  suitable  for  the  subsequmt 
analysis  of  micromotion  or  grain  rotation. 

The  preferred  system,  whether  cinema  or  video  cameras  are  used,  ^>pears  to  be  that  which 
uses  one  camera  per  specimen  window.  This  restriction  dictates  the  need  to  devise  a 
method  of  equating  a  given  image  with  loading  cap  displacement.  Towards  this  end,  a 
frame  counter  was  used  for  each  of  the  cinema  cameras  and  the  ai^ropriate  frame  number 
recorded  at  critical  points  in  the  loading  sequoice,  i.e.,  start  of  cap  displacement,  reversal 
of  cap  displacement  during  cyclic  loading,  and  maximum  displacemoit,  etc.  This  system  is 
suitable  only  if  the  loading  and  unloading  strain  rates  are  known  and  constant. 
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4.3. 1.2  Data  Capture  Through  Video  Recording.  An  alternate  technique  which  relies  on 
victoo  camera  imaging  was  investigated  for  quantifying  the  surficial  microstructural 
re^Kmse.  The  feasibility  of  using  such  "machine  vision”  to  quantify  microstructural  effects 
was  denumstrated  and  a  prototype  system  (hardware  and  software)  is  bdng  develq)ed. 
Preliminary  results  are  reported  to  aid  the  reader  in  choosing  the  appropriate  technique 
(either  cinematic  or  video  imaging)  for  subsequmt  research. 

The  developmental  video  system  uses  three  high  resolution  black  and  white  video  cameras,  a 
laser-disc  video  recorder,  and  a  sequencer  to  ensure  that  images  are  obtained  from  the 
sqipropriate  camera  at  the  proper  time.  The  principal  advantage  of  using  the  video  system  is 
the  ability  to  consistently  determine  the  centroid  and  orientation  of  a  given  grain  in  each 
frame.  Grain  fracture  and  disintegration  can  also  be  quantified  unambiguously  by  such  a 
system  given  suitable  image  analysis  software.  Because  these  microstructural  effects  are 
important  in  the  evolution  of  fabric  under  loading,  they  must  be  quantified  in  an  unbiased 
and  consistent  marmer. 

Several  processes  must  be  performed  on  a  single  digitized  image  of  the  iqjecimen  to  permit 
the  detection  and  tracking  of  individual  grains  and  remnants  from  image  to  image.  The 
images  must  undergo  preconditioning  to  maximize  the  contrast  between  adjacent  grains  and 
between  grains  and  voids.  The  high-contrast  image  must  then  be  "labeled”  to  sqnrate  the 
grains  which  are  being  tracked  from  the  remaindo’  of  the  image.  Subsequently,  grain 
centroids  and  orientation  must  be  determined  to  pmnit  comparison  between  their  current 
location/orientation  and  that  in  the  previous  image.  The  results  of  this  process  then  serve  as 
the  basis  for  analysis  which  includes  determination  of  average  direction  and  velocity  of 
motion  with  respect  to  the  fiducial  marks,  grain  rotation,  the  formation  of  clusters,  grain 
fragmentation,  and  the  growth  of  damage.  The  individual  tasks  in  this  process  are  detailed 
next. 

Preconditioning  is  performed  to  produce  the  best  possible  image  for  labeling.  The  three 
types  of  preconditioning  used  are  histogram  equalization,  median  filtering,  and 
homomorphic  Altering.  Histogram  equalization  is  a  commonly  used  contrast  stretching 
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technique  which  broadens  the  gr^  level  distribution  of  an  image  so  that  tte  image  sqians  the 
entire  grey  level  range  available.  This  technique  does  not  alter  the  information  in  the 
image.  Rather,  it  matches  the  image  to  the  range  of  grey  levels  best  perceived  by  humans. 
Hi^cgram  equalizaticm  can  be  aj^lied  to  the  image  as  seen  by  tlw  video  camera  before  a 
test  is  conducted.  This  aspect  oisures  that  rnily  high  contrast  images  will  be  provided  for 
analysis.  "Hot  spots*  due  to  uneven  lighting  can  be  minimized  by  this  technique. 

The  second  preconditioning  technique,  median  filtering,  is  a  process  used  to  reduce  noise  in 
the  image.  A  median  filter  is  a  nonlinear  spatial  filter  which  smoothes  out  raiulom  noise  but 
preserves  the  grain  edges. 

The  last  method,  homomorphic  filtering,  is  used  to  improve  the  evenness  of  the  image.  This 
filtering  technique  is  based  mi  the  multiplicative  relationship  between  reflectance  and 
illumination  which  constitutes  an  individual  pixel’s  brightness.  Brightness  equals  reflectance 
times  illumination.  Because  of  this  relationship,  reflectance  and  illumination  cannot  be 
filtered  separately  using  linear  filters.  Therefore,  one  must  take  the  logarithm  of  the  image, 
which  converts  pixel  brightness  to  the  form  log  (reflectance),  plus  log  (illumination).  One 
can  then  t.«ke  a  two^mensional  discrete  fast  fourier  transform  (2-D  FFT)  of  the  image  and 
tq)ply  a  frequency-domain  filter.  This  suppresses  low  frequencies  which  correspond  to 
general  illumination.  Taking  the  inverse  2-D  FFT  and  then  the  antilog  of  the  transform 
results  in  a  "reflectance-only"  image.  As  reflectance  is  a  material  prc^rty,  the  image  is 
necessarily  more  even. 

The  labeling  process  seeks  to  recognize  object  classes  in  the  image  (grains,  remnants,  voids, 
etc.)  based  on  spatial  and  grey  level  properties.  Tliis  process  results  in  similar  objects  in 
each  class  being  distinguished  by  a  unique  grey  level.  A  threshold  is  then  applied  to 
highlight  objects  of  interest,  either  grains,  remnants,  voids,  etc. 

A  morphological  filter  is  applied  following  preconditioning.  This  technique  selects  the 
objects  in  the  image  which  most  closely  conform  to  the  physical  features  dictated  to  the 
filter;  i.e.,  shape  and  size  of  the  grains,  typical  area  covered  by  a  grain  or  remnant,  etc. 
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Individual  gnin  and  du^  silhou^tes  have  been  ictendfied  at  this  point.  Following 
ayjfdication  of  the  mcvphdogical  filter,  determinatitm  of  the  grain  centroids  and  orientation 
is  a  rdadvely  simple  i»ooess.  As  indicated  eariier,  Uie  morphology  is  dictated  to  the  filler. 
This  teduiique  can  easily  be  adapted  to  recognize  a  wide  variety  of  geometric  slopes  or 
other  microstructural  aq)ects.  Thus,  this  automated  process  may  eventually  be  used  for  the 
microstructural/micromotion  analysis  of  a  wide  vari^  of  soils  such  as  well-graded  soils 
and/or  clayey  soils. 

Tracking  is  the  process  by  which  the  centroids  and  orientation  of  selected  grains  are 
dmected  and  recorded  firom  image  to  image.  Equally  important  is  the  added  cspability  of 
monitoring  the  growth  of  damage  within  the  viewed  field  through  the  implementation  of  an 
adaptive  moipholt^cal  filter.  Such  a  filter  will  recognize  the  existence  of  remnants, 
may  have  dramatically  different  geometric  pnperties  than  the  parent  grain,  and  continue  to 
track  their  micromotions/oiientations. 


Analysis  involves  compilatitm  of  the  infbrmaticm  obtained  during  tracking  to  highlight  a 
desired  feature  of  the  microstructural  response.  Clustering  and  rfamagf.  growth  are  two 
important  features  which  are  being  examiiwd.  The  consistency  and  precisimi  associated  widi 
the  automated  processing  and  analysis  of  images  provide  a  d^ree  of  accuracy  whidi  cannot 
be  efficiently  attained  through  manual  analysis. 


4.3.2 


Recording  Using  X-Rav  Techniques 


Several  tests  were  conducted  using  a  120-keV  x-ray  device  (Ref.  36).  These  tests  involved 
the  placement  of  1.3-mm-diam.  lead  pellets  (called  particles)  within  a  specimen  of  O.Sfi-mm- 
diam.  sOica  grains  in  the  array  shown  in  Figure  21.  The  specimen  container  was  configured 
to  have  two  opposing  windows  at  the  same  height  (0.56  h),  thus  permitting  transmission  of 
the  X  rays  through  the  specimen.  Three  tests  were  conducted  with  macroscopically  similar 
specimens  being  subjected  to  6.0-percent  strain.  The  test  r^me  (macroscopic  strain 
application  and  removal  of  0.5  potent  per  minute)  was  similar  to  that  used  during  tests 
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a.  Loading  csqp  di^laceiTMmt  (mm). 


b.  Particle  anangement.  c.  Zone  definititHi. 

Figure  21.  Displacement  of  avoage  particle  embedded  within  upper/ 
miMe/lower  zones  of  viewing  window  (window  at  0.S6  h). 
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in  «4iic^  dnemt  imnges  were  obtained.  These  tests  provided  qualitative  insight  inio  the 
influence  oi  sidewall  flktkxi  on  the  surfldal  grain  motioo  conqared  to  the  motion  of 
particles  located  in  the  center  the  qiecimen  which  are  presumably  unafie^  by  sidewall 
friction.  Although  the  use  of  x  nys  cannot  yet  provide  the  detailed  images  oi  paitide 
motioos  (displacement  and  rotation)  which  are  seen  in  the  cmema  images,  conqmison  oi  die 
‘average  grain*  motion  (determined  by  tracking  die  motioos  ci  10s  of  grains  widiin  a 
predetermined  zone)  provides  a  means  of  inferring  whmher  internal  nuCTomotioos  are  the 
same  as  die  surfldal  mkromotions. 


4.3.3 


lOtBQIgL 


il  Response 


The  internal  microstructural  leqxxise  is  the  changes  which  occur  within  the  spedmen  as 
observed  at  discrete  stages  during  the  i^lication  of  macroaco|nc  uniaxial  strain.  This 
response  is  limited  to  changes  ^di  can  be  dMerved  in  die  qiedmen  after  it  has  been 
loaded  to  some  incranent  of  macroscopic  strain  (2  percent,  4  percent,  6  percent,  elc.)  and 
then  unloaded.  Miermnotions  cannot  be  observed  by  examining  qiedmens  which  are  not 
subject  to  a  changing  stress  (strain)  state  during  observation.  Rather,  observadoos  are 
limited  to  a  study  of  changes  in  grain  mentation,  damage  levels,  and  the  distribution  of 
damage  diroughout  the  ^ledmen.  Inftmnation  regarding  diese  changes  in  the  ^ledmen  was 
obtained  using  an  epoxy  impr^nadon  tedmique  with  subsequent  sectioning  to  expose  the 
interior  of  die  qiedmen.  These  techniques  are  (fescribed  next 

4.3.3. 1  Eooxv  Imorepiadon.  Wafering.  and  Polishing.  Standard  petrographic  procedures 
were  used  for  the  epoxy  impregnadon  of  the  tested  ^ledmens  as  wdl  as  the  subsequent 
sectioning  and  polishing  of  the  exposed  cut.  Aiehlo’,  Ltd.,  Lake  Bluff,  Illimns,  offers  a 
variety  of  systems  which  are  suitable  for  this  type  of  microstructural  analysis.  As  reladvely 
large  ^ledmens  were  being  examined,  rdined  techniques  had  to  be  developed  fev  eadi 
phase  of  the  microstructural  invesdgadon.  Widi  r^ard  to  impr^nadon,  a  low  viscosity 
resin  was  used  as  the  epoxy  mixture  had  to  flow  eadly  within  the  voids  betweoi  the  grains. 
The  mixture  had  to  be  "worked"  into  the  specimen  by  subjecting  the  specimen  to  a  vacuum 
while  it  was  submerged  in  an  epoxy  bath.  The  vacuum  served  to  draw  air  out  of  the  vdds 
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in  tile  speramen  which  wm  then  filled  with  epoxy  mixture  when  the  specimen  was  acain 
ntige^  to  atmoqiheric  pressure.  Several  cycles  of  "vacuumiog”  were  necessary  to 
adueve  full  inquegnsdion.  As  tiie  duration  and  number  of  vacuum  cycles  are  dqiendent  on 
the  void  siae  and  vdume,  epoxy  tenqierature  and  viscosity,  and  qiecimen  sire,  the 
impregnation  r^ime  is  best  dmremined  through  trial  and  error.  Extra  care  must  be  taken  to 
keq>  tile  specimen  and  qwxy  at  tiie  low  end  the  woridng  temperature  fi»  tiie  chosen 
qpoxy  because  rd^vely  lugh  temperatures  cause  the  epoxy  to  cure  more  iiqiidly.  Rapid 
curing  results  in  eitiier  shrinkage  cracks,  as  a  minimum,  or  the  onset  of  an  exothermic 
reactioi  witiiin  the  qiedmen  which  causes  uncontrolled  curing.  The  qwdmen  is  destroyed 
in  eitiier  case  as  the  fabric  is  distorted  and  pievioiisly  intact  (possibly  fractured)  grains  are 
broken  into  many  pieces.  Damage  occurs  due  to  impregnation  indiich  might  wrongly  be 
attributed  to  damage  whidi  was  supposed  to  hare  occurred  during  loading. 

The  Sfi-mm-sire  (cubic)  qiedmen  ippeared  to  be  tiie  largest  qiedmen  ndiich  could  be 
sectioned  by  ccmimetciaUy  available  diamond  saws  (assuming  the  saw  is  not  modified  to 
accommodate  larger  specimens).  Following  inqxrq;nation,  the  specimen  must  undergo 
additional  preparation  to  remove  excess  epoxy  from  the  sides  and  top.  This  is  accomidished 
by  either  cutting  or  grinding  and  serves  to  restrxe  tiie  qiedmen  to  its  nominal  dimensions  of 
50  mm  on  a  side.  The  specimen  can  tiien  be  halved  and  subsequently  sectioned  into  slices 
and  smaller  pieces  as  appropriate. 

R^ardlres  of  the  sire  and  orientation  of  the  exposed  face,  grinding  and  polishing  were 
necessary  to  prqiare  the  face  for  tptical  analysis.  Several  commercially  available  systems 
are  capable  of  producing  a  surface  suitable  for  tptical  analysis,  thus  techniques  for  use  witii 
a  pedfic  system  will  not  be  discussed  here.  The  difficulty  of  obtaining  a  highly  pdished 
surface  is  closely  related  to  the  amount  of  damage  the  pecimen  has  undergone.  This  is 
mainly  because  the  epoxy  may  not  fully  impregnate  grains  which  have  fractured  but  not 
shattered.  These  grains  tend  to  break  and  crumble  as  polishing  inogtesses.  The  remnants 
are  subsequently  drawn  across  the  face  of  the  pecimen  and  scarring  results.  Such  scarring 
may  interfere  with  the  ability  of  most  off-the-dielf  optical  analysis  systems  to  discern 
mioostructural  effects  of  intmest;  e.g.,  damage  and  fracture  of  individual  grains.  It  is  often 
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a  trade-off  between  the  amount  of  inftmnation  <me  wishes  to  obtain  from  an  image  and  the 
ability  of  the  optical  systmn  to  perform  the  analysis  as  fine  polishing  is  required  to  eaqxne 
smne  of  the  more  subtle  damage  occurring  in  the  spedmea.  More  importantly,  the  visibte 
(bunage  may  be  more  a  functkm  of  the  polishing  jffocess  than  a  function  of  die  true  damage 
existing  in  the  qiecimen.  One  must  also  consider  the  fact  that  the  peak  load  has  been 
removed  and  the  specimen  allowed  to  relax  before  impr^nadon.  Apfnoximately  30  percent 
elastic  rebound  of  the  local  disidaoement  occurs  during  load  removal.  Thus,  the  grains  are 
collectively  under  a  relatively  low  stress  at  die  time  of  grinding  and  polishing,  yet  ^laliing 
still  occurs.  This  is  worth  noting  because  of  the  expected  difficulty  in  obtaining  a  well- 
polished  fiice  on  a  qwdmen  imiv^nated  while  still  under  load.  Although  such  an  image  is 
desirable,  the  amount  of  spalling  and  scarring  might  preclude  obtaining  a  useful  image. 
Techniques  have  been  amceived  to  impr^nate  ^)ecimens  un^  load;  however,  they  have 
not  been  pursued. 

4.3.3.2  Optical  Analysis.  An  Olympus  CUE  system  was  used  for  optical  analysis 
purposes.  Several  ciqnbilities  exist  within  the  system  to  aid  in  the  dderminatitm  of 
microstructural  dianges  within  the  ^ledmen.  These  algorithms  include  edge  detection, 
small  feature  enhancement,  various  filters  for  highlighting  the  boundaries  of  grains,  and  the 
ability  to  quantify  a  limited  number  of  geometric  aq)ects  of  individual  grains  and  their 
(nientatitm.  The  suitability  of  these  c^abilities  and  their  use  in  the  development  of  a 
microstructural-based  omstitutive  will  be  discussed  in  Subsection  6.4. 
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S.O  SPECIMEN  PREPARATION  AND  TESTING 

S.l  EBBEARAUQN 

All  tests  woe  conducted  using  spectmens  manuCKtuied  from  either  O.SO-mm-  or  0.7S>mm' 
diam.  silica  grains  (7  »  2.65).  Macroscopically,  each  qwcinien  had  a  dry  unit  weight  of 
1.7  g/cc.  Thus,  the  v<^  ratio  (e)  was  0.61  percent  and  the  porosity  (n)  was  38  percent. 

A  measured  weight  of  grains  was  unifrnrmly  rained  into  the  qieciinen  container  and 
subsequently  vibrated  until  the  desired  maatrscopic  propeities  were  obtained.  Hie  rigid 
loading  cap  (Subsection  4.1)  was  then  jdaced  atop  die  soil  and  die  assembly  jdaoed  within 
die  loading  frame  for  testing. 


All  tests  were  diaplaoement-cewitrolled.  The  diiyilaceinent  rate  equated  to  a  macroscopic 
loading/unloading  rate  of  O.S>percent  strain  per  minute.  The  two  types  of  tests  conducted 
were  monotonic  loading  and  unloading  as  wdl  as  cyclic  loading  with  hi^ier  macroscopte 
strain  being  applied  with  eroh  cycle.  Unloading  was  typically  initiated  immediately  upon 
attaining  the  desired  macroscopic  strain,  but  the  peak  strain  (peak  displacement  of  the 
loading  cap)  was  hdd  constant  after  some  peak  macroscopic  strain  aiplicatimis  to  investigate 
the  stress  relief  properties  of  the  granular  material.  The  test  regime  is  presented  in  Table  2. 
The  first  two  digits  in  the  test  number  designate  the  grain  size  (SO  0.50  mm, 

75  «  0.75  mm),  the  secemd  two  designate  die  peak  macrosctpic  strain,  and  the  characta* 
within  the  parentheses  derignates  whether  strain  application  was  monotonic  (M)  or 
incremrotal.  If  the  strain  was  incrementally  rpplied,  the  number  in  the  parentheses  indicates 
the  size  of  the  increment.  Test  designations  with  a  suffix  of  I  or  n  indicate  that  like  tests 
were  performed  with  srpphire  windows  at  different  heights  to  determine  micromotimis  at 
varkHis  locadmis  on  the  side  of  tb  ximen. 
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Table  2.  Uniaxial  strain  ai^lication  test  regime. 


Test 

Macroscc^ic 
Strain  Regime 

Window 

Height  (h) 

Peak 

Displacement 

Retained 

50.03  (M) 

0-3-0 

0.56 

No 

50.06  (M) 

0-6-0 

0.56 

Yes;  6.5  min 

50.09  (M) 

0-9-0 

0.56 

Yes;  7.5  min 

50.09  (3) 

0-3-0-6-0-9-0 

0.56 

Yes;  1.0  min  at 
each  peak 

75.05  (M) 

0-5-0 

0.375 

0.500 

0.750 

Yes;  4.0  min 

0-1-0-2-0-3- 

0-4-0-5-0 

0.750 

Yes;  4.0  min 
at  5.0%  strain 

75.06  (M) 

0-6-0 

0.56 

No  (Interior 
particle 
displacement 
measured  using 

X  rays.) 

75.10  (M)  I 

0-10-0 

0.375 

0.500 

Yes;  5.0  min  at 

10%  strain 

75.10  (M)n 

0-10-0 

0.500 

0.750 

Yes;  5.0  min  at 

10%  strain 

75.10  (2)1 

0-2-0-4-0-6- 

0-8-0-10-0 

None 

Yes;  4.0  min  at 

10%  strain 

75.10  (2)  n 

0-2-0-4-0-6- 

0-8-0-10^ 

Yes;  4.0  min  at 

10%  strain 

NOTE:  Suitable  photographic  records  of  grain  displacement  and  fnu:ture  were  only  obtained 
for  0.75-mm-diam.  grains  as  the  O.SO-mm-diam.  grains  proved  too  small  to  be 
manually  (and  precisely)  distinguishable  from  frame  to  frame.  Therefore,  only  the 
analysis  results  of  the  0.7S-mm  grain  diameter  cinema  records  will  be  present^. 
General  trends  observed  in  the  micromotions  of  the  O.SO-mm-diam. 
grains  are  included  in  the  interpretation  of  the  results. 
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6.0  TEST  RESULTS  AND  DISCUSSION 


Two  ^jpfoaches  were  used  in  evaluating  the  microstructund  effects  which  woe  induced  and 
observed  during  the  s^lication  of  uniaxial  ^rain.  Tte  rirst  ^jfHoach  invcdved  reccnding  die 
micromotions  and  observation  of  other  microstructural  effects  which  occurred  on  the 
qiecimen  surface.  Such  effects  included  fracnire  and  disintegration  of  grains  and  the 
formation  of  damage  chains  and  clusters.  The  seccmd  tqiproach  involved  q»xy 
impregnatitm  of  unloaded  specimens  which  were  dien  secdtmed  and  polished  to  reveal 
changes  in  the  internal  microstructure  of  the  specimen.  The  results  and  discussion 
associated  with  each  aj^roach  will  be  presented  in  sqiarate  subsections.  Subsectitm  6.2 
presents  the  test  results  in  macroscopic  terms,  i.e.,  uniaxial  stress-strain  reqxmse,  pretest 
and  posttest  grain  size  distribution  curves,  etc.  Subsection  6.3  presoits  the  surficial 
microstructural  refuse  and  Subsecticm  6.4  presents  the  internal  microstructural  reqxmse. 


W :  Vi  a  I I L* J  J ft  I ;  1 4^1  g  1 3 1  a C ■ W  C«i  it  J m : 


Table  3  summarizes  the  peak  axial  stress  associated  with  peak-apfdied  (incremental  or  final) 
strain  for  each  test.  Figures  22-30  present  the  stress-strain  curves  associated  with  each  test 
and  the  development  of  axial  stress  versus  time. 

Figures  31-33  show  the  grain  size  distributitMi  charts  for  virgin  materials  used  in  the  t^ 
series  as  wdl  as  for  those  materials  subjected  to  the  same  macroscqiic  strain  levels 
(incremental  and  peak)  as  the  tested  specimens.  The  first  two  digits  in  the  grain  size 
distribution  curves  indicate  the  maximum  macrostxrpic  strain  which  was  tyrplied  to  the 
qredmen.  The  character  within  the  parentheses  designates  whether  strain  applicatim  was 
M  or  incremratally  applied.  If  the  strain  was  incrementally  allied,  the  number  in  the 
paraitheses  indicates  the  size  of  the  increment.  For  example,  a  test  designated  06(02) 
indicates  that  a  specimen  was  subjected  to  a  peak  macroscopic  strain  of  6.0  percent  in 
increments  of  2.0  percoit.  As  always,  the  macrostrain  was  removed  from  the  specimen 
before  the  application  of  the  next  strain  increment. 
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e  3.  Summary  of  axial  stress  fw  given  incremental  or  peak  strain. 
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NOTE:  ‘Axial  stress  is  given  in  units  of  megapascal  and  kips  per  square  inch. 


STRFSS  (MPA)  AXIAL  STRESS  (MPA) 


STRESS  (MPA)  AXIAL  STRESS  (MPA) 
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eradiation  curves  for  0.75<inm-dUun.  particle  specimnu  (virgin,  04[2]  and  10(2]) 
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Figure  32.  Gradiation  curves  for  0.7S-imn-diam.  particle  q)ecimens  (virgin,  03IMJ,  03[11). 
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Figure  33.  Gradiation  curves  for  O.SO-iran-diam.  pvticle  q)eciiiiens  (virgin,  03[M],  06[M],  09[M].  0913]). 
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The  listed  results  indicate  nothing  unusual  in  the  expected  response.  Generally  both  types  of 
qweimens  (SO-  and  7S-mm  grains)  denKxistiated  a  rdatively  weak  Young’s  modulus  during 
initial  strain  tq^lication  (0  to  1.0  pocent  macrosct^c  strain)  with  the  strength  of  the 
material  increasing  as  strain  was  applied.  This  stiffening  occurred  from  1.0-  to 
~3.0-pacent  macroscopic  strain  for  the  O.SO-mm-diam.  grain  q)ecimens  and  ftom  1.0-  to 
-S.O-percent  macroscopic  strain  for  the  0.7S-mm-diam.  grain  ^)ecimens.  It  appeared  dtat 
after  »5.0  percent  macroscqnc  strain  had  bera  sqjjplied  to  a  givra  specimen,  a  slightly 
decreased,  yet  cemstant.  Young’s  modulus  was  exhibited. 

6.3  SURnClAL  MICRQSIRUCTURAL.  RESgQNSE 

As  mentiemed  previously,  the  surfleial  microstructural  req)onse  considm  the  micromotions, 
fracture,  and  disintegration  of  individual  grains  tu  predetermined  locations  on  the  surface  of 
the  ^)ecimai.  The  following  observations  were  made  from  cinema  records  obtained  through 
the  sapphire  windows.  Two  approaches  were  used  in  the  interpretation  of  the  images. 

Pint,  the  motion  of  the  "average  grain"  within  the  viewed  neighborhood  was  considoed. 
Secemd,  the  images  were  viewed  with  the  intent  of  develq>ing  an  understanding  of  the 
growth  of  clusters,  the  motion  of  such  clusters,  and  the  develc^ment  of  damage  within  the 
specimen.  Each  of  these  responses  will  be  discussed  because  of  the  limitation  imposed  by 
the  presence  of  the  sapphire  viewing  window.  A  qualitative  description  of  the  grain 
micromotions  and  the  evolution  of  damage  is  givoi  next. 

The  following  observations  were  made  through  the  window  at  0.56  h  with  specimois 
entirely  made  up  of  either  0.50-  or  0.75-mm-diam.  grains.  It  appears  that  individual  grains 
shift  and  reorient  themselves  during  the  early  stage  of  strain  2q>plication  (<  1.0  percent 
macroscopic  strain)  and  remain  intact.  As  strain  application  progresses,  the  stress  at  grain 
contacts  for  the  now  "locked  in"  grains  increases  until  the  stress  level  is  adequate  to  fracture 
and  crush  individual  grains.  Such  stresses  appear  to  have  fiiUy  developed  by  the  time  the 
0.50-mm-diam.  grain  specimens  have  been  subjected  to  3.0-percent  macroscopic  strain  and 
the  0.75-mm-diam.  grain  specimens  have  been  subjected  to  5.0-percent  macroscopic  strain. 
After  these  levels  of  macroscopic  strain  have  been  attained,  fracture  and  crushing  of 
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indivklual  grains  are  readily  apparent.  Even  though  t)»se  observations  are  made  from 
viewmg  trough  the  window  at  0.56  h,  the  same  general  obsm'vations  woe  macte  dirough 
windows  at  other  heights.  It  is  reasonable  to  assume,  therefore,  that  each  of  these  reqxmses 
coexist  in  the  qiecimai  during  any  stage  of  load  tqpplkation.  Even  at  low  macroscopic 
strain  levels,  combinations  of  these  factors  may  result  in  highly  tocalized  damage.  Mt»e 
importantly,  it  is  reasonable  to  assume  that  all  of  these  responses  are  occurring 
simultaneously  in  the  specimen  during  loading.  These  responses  continue  to  occur  in 
spedmois  held  at  some  peak  level  of  macroscopic  strain  as  the  resisting  load  decreases  with 
time  after  the  peak  strain  level  has  been  attained.  The  mechanisms  of  individual  grain 
displacement,  rotation,  fracture,  and  di^tegration  can  thus  be  considered  the  principal 
causes  for  time-dq)endent  strain  (creep)  under  constant  load. 

With  regard  to  the  modem  of  the  "average  grain,"  die  individual  vertical  and  horizontal 
displacements  of  selected  grains  within  a  given  cinema  frame  (approximately  SO  grains  per 
series  of  frames)  were  summed  with  the  resulting  average  motion  being  assigned  to  an 
imaginary  average  grain.  The  motions  of  these  average  grains  are  shown  in  Figures  34-37. 
Although  one  would  expect  the  averaged  grain  motion  to  appear  rather  continuous  and  be  a 
smooth  curve,  the  path  is  erratic  and  nonlinear.  Furthermore,  the  path  of  the  average  grain 
does  not  retrace  the  path  it  followed  during  the  loading  phase.  The  wide  variation  in  this 
re^nse  is  evidoit  in  Figure  36,  which  shows  the  avoage  grain  rdxHmding  more  than 
expected  in  one  case  (viewing  window  at  0.75  h)  and  less  in  another  (viewing  window  at 
0.50  h)  within  the  same  specimen.  In  this  test,  the  unloading  paths  of  the  average  grains 
veer  toward  the  linear  displacemoit  path  suggested  by  a  local  theory  formulation.  Two  facts 
are  evident.  First,  although  one  might  expect  the  loading  and  unloading  paths  to  be 
nonlinear  and  lag  behind  the  predicted  linear  displacement  path  due  to  wall  friction,  this  is 
not  the  case.  Second,  such  microstructural  responses  remain  hidden  if  one  considers  only 
the  macroscopic  material  response.  Regarding  the  loading  and  unloading  paths,  wall  friction 
is  considered  to  retard  the  vertical  motion  of  individual  and  groups  of  grains.  If  this 
influence  was  considered  significant,  the  motion  of  all  grains  adjacent  to  the  walls  or 
windows  would  lag  behind  the  displacements  predicted  by  assuming  a  linear  displacemoit 


73 


PL-TR-92-1064 


Macroscopic  Vertical  Strain 


Figure  35.  Displacement  of  average  particle  for  observations  made  at  0.75  h; 
cyclic  loading  to  1.0-,  2.0-,  and  3.0-percent  macroscopic  strain, 
monotonic  unloading  between  cycles;  and  0.75-mm-diam.  particles. 
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Figure  36.  Displacement  of  average  particle  for  observations  made  at  0.375, 
0.50,  and  0.75  h;  monotonic  loading  to  6.0-percent  macroscopic 
strain,  monotonic  unloading;  and  0.75-mm-diam.  particles. 
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gradient  throughout  the  specimen.  Indeed,  di^lacement  at  some  windows  exceeded  those 
predicted  by  a  local  di^lacement  theory.  One  can  conclude  that  the  sidewall  friction  does 
influence  the  micromotion  of  grains  adjacent  to  the  wall;  howeva-,  the  influence  is  not  so 
great  that  the  natural  mobility  of  the  grains  is  destroyed.  It  is  2q)parent  that  the 
microstructural  responses  may  provide  variatitH)  around  the  "theme”  presented  by  observing 
the  macrostructural  resptmse.  This  variation  manifests  itself  in  regions  of  relativdy  low  and 
high  displacement,  strain,  grain  rotation  and  damage,  etc.  All  of  these  responses  occur  in 
various  measures  during  the  loading  process  and  contribute  to  the  macrostructural  response. 
The  evolution  of  and  interplay  between  these  reponses  are  obviously  complex  and  are  most 
likely  the  basis  for  variations  in  the  results  of  tests  using  macroscopically  similar  specimens. 

Some  measure  of  the  variation  in  the  vertical  of  grains  in  a  specimen  subjected  to  uniaxial 
strain  is  given  in  Figure  38.  This  figure  shows  a  summary  of  the  vertical  diplacements  of 
some  650  grains  as  observed  through  a  window  at  0.S6  h  with  the  0.7S-mm-diam.  grain 
specimen  while  subjected  to  3.0-percent  macroscopic  strain.  The  expected  vertical 
diplacement  of  the  grains  (if  a  linear  diplacemoit  gradient  is  assumed)  is  shown  by  the 
v^cal  line  designated  local  theory  prediction.  However,  some  of  the  grains  are  immobile 
and  some  have  diplaced  pproximately  three  times  the  expected  distance.  It  is  important  to 
note  that  the  average  measured  diplacement  (each  grain  having  an  equal  "Weight”)  is  equal 
to  the  diplacement  predicted  by  the  assumption  of  a  linear  diplacement  gradient.  This 
variation  in  the  mobility  of  grains  within  a  grouping  or  "neighbortiood"  provides  insight  into 
the  dissipation  of  energy  pread  by  the  specimen  as  macroscopic  strain  is  applied. 
Additionally,  this  variation  serves  as  the  basis  for  the  development  of  a  theory  which  uses 
the  diffusion  equation  for  the  prediction  of  micromotions  within  a  specimen  subjected  to 
uniaxial  strain.  This  theory  is  presented  in  Reference  40. 

The  relationship  between  vertical  and  horizontal  grain  diplacements  as  a  function  of 
macroscopic  strain  is  shown  in  Figure  39.  The  micromotions  were  viewed  through  a 
window  at  0.S6  h.  Data  for  this  figure  were  produced  by  summing  the  vertical  motions  of 
all  grains  within  a  series  of  frames  at  predetermined  increments  of  macroscopic  strain. 
Likewise,  the  horizontal  motions  of  all  grains  were  summed  over  the  same  series  of  frames. 
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Figure  38.  Variation  in  the  vertical  displacement  of  sand  grains  under  uniaxial  strain 

(0.75-mm-diam.  particles,  window  at  0.56  h,  3.0-perccnt  macroscopic  strain). 


Horizontal  Displacemont 
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Reli.  ive  Horizontal  Displacement  (Z) 
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Percentage  of  vertical  and  horizontal  di^lacement  versus  an>lied 
macroscopic  strain. 
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The  di^byed  in  Figure  39  shows  the  relative  amounts  of  average  vertical  and 
horizontal  dispbcements  with  reflect  to  the  total  observed  di^lacement.  A  cyclical 
relatkmship  betwera  the  magnitudes  of  the  vertical  and  horizontal  di^laoements  aj^iears  to 
devdop  as  uniaxial  strain  is  applied.  The  grains  move  primarily  vertically  downward  during 
the  tq)plicati<m  of  between  0-  and  2>percent  macroso^ic  strain.  As  strain  plication 
continues,  the  grains  move  laterally  into  a  more  compact  arrangement.  This  lateral  motkm 
occurs  between  2-  and  4-percent  macroscopic  strain.  At  that  point,  grain  fracture  becomes 
significant.  The  cycle  then  rq)eats  itself  as  the  grains  and  remnants  di^lace  vertically  and 
then  shift  horizcmtally  into  an  even  more  compact  arrangement  before  fracturing  again. 

Figure  40  shows  the  tracks  made  by  grains  which  di^lace  along  the  sidewalls  of  the 
confinement  vessel.  The  photogrq)hs  were  taken  following  the  application  of  lO-pocent 
macroscopic  strain.  Grains  and  their  remnants  tend  to  move  towards  the  nearest  bottom 
ccrnier  of  the  specimen  container.  Lateral  motion  of  the  grains  in  the  lower  half  of  the 
q)ecimen  is  most  pronounced  with  grains  only  alcmg  the  wall  center  liiw  exhibiting  generally 
vertical  displacemoits. 

The  growth  of  clusters  and  development  of  damage  are  two  tesptmses  which  must  be 
considered  jointly  as  one  examines  the  surfidal  microstructural  reqxmse.  As 
previously,  the  grains  displace  vertically  in  response  to  the  sy>plication  of  macroscopic 
uniaxial  strain.  This  resptmse  predominates  until  the  stress  in  individual  grains  becomes 
significant  enough  to  result  in  fracture  of  the  grain.  The  fracture  of  these  grains  iqjpears  to 
occur  randomly  during  the  onset  of  damage.  Ho^i^er,  as  additional  macrosct^c  strain  is 
applied,  chains  of  fractured  grains  appear  to  grow  b^een  the  grains  which  fracture  first. 
Such  chains  are  highly  evident  if  the  specimen  is  subjected  to  >  7.0-percent  macroscofnc 
strain.  These  damage  chains  are  considered  to  be  the  boundaries  of  clusters  within  whidi 
there  is  relatively  little  damage.  Clusters  are  three-dimensional  structures  which  are  formed 
within  the  specimen.  Thus,  only  a  portion  or  section  of  a  cluster  is  visible  through  the 
ssqqrfiire  window.  Full  definition  of  the  cluster  geometry  and  its  evolution  unctor  uniaxial 
strain  dictates  the  ability  to  "see”  into  the  specimm  as  it  is  loaded.  This  is  impossible  using 
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Figure  40.  Grain  displacement  Uacks  on  the  inside  of  the  confinement  vessel  walls. 
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Figure  40.  Continued. 
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(d)  Specimen  Mde  of  loading  ptate. 

Figure  40.  Concluded. 
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the  cuncot  experimental  setup.  Furthennore,  the  micromotion  of  a  well-defined  cluster  is 
impossible  to  discmn  usi  .g  the  current  s  "Jp.  It  is  thought  that  at  macroscopic  strain  >5.0 
percent,  the  micromotioi  and  intmactiow.>  of  the  clusters  may  influmux  the  specimen 
macrcMcopic  response  more  than  interactions  between  individual  grains.  Additionally,  die 
geometry  of  the  clusters  may  be  such  that  they  efficiendy  mobilize  the  stroigth  needed  to 
resist  the  continued  t^if^ication  of  uniaxial  strain.  Subsequent  rotadmi  of  the  principal  stress 
axes  following  the  triplication  of  the  initial  increment  of  uniaxial  strain  will  typically  roida* 
a  stress-strain  response  significandy  different  than  die  response  of  a  virgin  qiecimen. 

Again,  investigation  of  this  .esponse  is  not  possible  using  die  current  device.  This  is 
because  one  is  limited  to  actively  viewing  only  the  respmise  of  surficial  grains,  and  the 
confinement  vessel  does  not  permit  rotatimi  of  the  principal  stress  axes. 

The  last  point  to  be  addressed  is  the  late-time  microstnictural  response  of  a  ^ledmen  which 
has  been  subjected  to  some  increment  of  uniaxial  strain  whidi  was  then  held  constant.  The 
peak  strength  (Figs.  23-30)  mobilized  by  the  qiedmen  immediately  decreases  upon 
attainment  of  the  peak  strain  level.  The  peak  mobilized  strength  decreases  because^ 
microstructural  activity  continues  after  the  peak  macroscopic  strain  level  has  been  attained. 

It  is  surmised  that  at  the  instant  that  the  final  increment  of  macroscopic  strain  is  aqiplied, 
there  is  some  finite  number  of  grains  which  are  on  the  verge  of  fracture.  Some  do  fracture 
and  some  may  not.  Adjacent  grains  receiving  Hot  load  may  or  may  not  fracture.  This 
process  continues  until  the  relatively  high  energy  concoitrated  within  several  ”key”  grains 
ends  up  being  more  evenly  distributed  throughout  the  qiedmen.  There  is  some  levd  of 
distribution  at  which  the  grains  and  their  remnants  can  withstand  the  allied  strain  with  no 
additional  fracture.  This  is  a  point  of  equilibrium.  This  premise  was  explored  by  analysis 
of  micromotions  following  attainment  of  the  peak  macroscopic  strain;  however,  little  late 
time  activity  was  recorded.  It  is  further  surmised  that  either  very  few  key  grains  fracture 
during  the  process  of  distributing  stored  energy  or  the  grains  which  fracture  reside  inside  the 
specimen.  The  rate  at  which  the  specimen  can  distribute  this  energy  may  be  an  important 
consideration  in  the  rate  effects  demonstrated  by  the  mechanical  response  of  some  soils. 
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The  micfomotion  of  the  surficial  grains  is  closely  tied  to  the  internal  microstructural 
raponse.  Prq)er  intnpretation  of  the  surficial  response  requires  additional  information  with 
r^ard  to  the  internal  inicrostrw:tUTal  respcmse  of  the  grain  assemblage.  The  requited 
informatitm  is  provided  by  x-ray  tracking  of  lead  particle  di^laceimnt  and  epoxy 
impr^naticm  and  sectitming  of  the  specimen. 


6.4 


ICii 


5RNAL  MICROSTRUCTURAL  RESULTS 


The  two  techniques  used  to  explore  the  internal  microstructural  respmse  of  the  ^)ecimens 
were  x-ray  tracking  of  lead  particles  and  epoxy  impr^nation  and  sectioning.  The  x-ray 
tracking  technique  is  used  during  the  application  of  macroscopic  strain  and  provides  the 
basis  for  comparison  of  surficial  and  internal  micromotions.  The  epoxy  impregnatimi 
technique  permits  the  direct  observation  of  microstructure  after  the  load  has  been  removed 
from  the  specimen. 

The  displacements  of  embedded  lead  particles  were  compared  to  those  of  the  surficial  grain 
displacements  in  Figure  41.  (The  diameter  of  the  embedded  lead  particles  was  1.3  mm 
whereas  the  nominal  sand  grain  diameter  was  O.SO  mm).  The  given  displacement  history 
shows  the  average  displacement  of  20  embedded  particles.  As  observed  through  a  window 
at  0.56  h,  posttest,  the  lead  particles  had  indentations  at  each  of  the  contacts  with  the  silica 
grains.  Each  lead  particle  had  10-12  indoitadons  whereas  the  average  coordination  number 
of  the  silica  grains  appeared  to  be  between  3  and  5.  As  a  result,  the  lead  particles  were 
more  rigidly  confined  than  the  silica  grains.  Additionally,  the  lead  particles  did  not 
fracture.  The  lead  particles  were  constrained  from  displacing  in  discrete  "jumps”  whereas 
the  silica  grains  were  not.  Thus,  their  displacement  paths  appeared  more  continuous.  The 
average  displacement  path  of  the  embedded  lead  particles  indicates  the  same  nonlinear 
di^lacemmt  path  during  both  loading  and  unloading  as  shown  by  the  average  grains  at  the 
specimen  surface. 
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Figure  41.  Vertical  displacement  of  average  particle  centroid  in  response  to 
^plied  macroscopic  strain. 
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However,  dqiaiture  of  the  displacement  path  from  that  predicted  by  local  theoiy  is  much 
less  pronounced  for  the  embedded  lead  particles  than  for  the  grains  on  the  specimen  surface. 
Thus,  there  appears  to  be  some  variation  between  the  displacement  paths  of  particles  in  the 
cento*  of  the  specimen  as  compared  to  those  which  are  adjacoit  to  the  sapphire  window. 

The  variatitm  is  found  mainly  in  the  magnitude  of  displacements  but  not  the  troids.  No 
informatimi  was  obtained  r^arding  the  distributitm  of  the  micromotitms  due  to  the  small 
number  of  embedded  lead  particles  (20).  It  appears,  though,  that  the  motion  of  the 
embedded  particles  is  more  of  a  continuous  nature  than  the  motion  of  grains  located  at  the 
boundary  of  the  specimen.  This  apparent  o}ntinuity  of  displacement  for  embedded  particles 
may  result  from  (1)  the  particles  are  free  to  displace  in  any  direction  whereas  the  surfoce 
grains  are  restrained  by  the  presence  of  the  ctmfinement  vessel  walls,  and  (2)  the  embedded 
particles  are  made  of  lead  and  are  larger  than  the  sand  grains.  The  embedded  particles, 
having  the  freedom  to  displace  laterally  in  any  direction,  have  inho-ently  greater  mobility 
than  the  grains  restrained  by  the  contains  walls.  It  is  the  restraint  which  leads  to  grain 
fracture  and  disint^ration  seen  through  the  sapphire  windows.  The  embedded  particles 
displace  and  rotate  more  subtly  than  the  wall-restrained  grains. 

While  grain  fracture  and  disint^ration  may  be  readily  apparoit  (hi  the  specimen  boundary, 
interior  grains  may  remain  relatively  intact  and  yet  retain  the  freedom  of  modem  to  orient 
themsdves  into  a  denser  arrangement  (inior  to  inevitable  fracture  and  disint^ration).  The 
foct  that  the  embedded  particles  are  made  of  lead  and  are  larger  than  the  sand  grains  also 
makes  their  average  motion  differ  from  that  of  the  sand  grains.  Being  made  of  a  malleable 
material,  the  lead  particles  simply  deform  in  reqmnse  to  the  forces  rqpplied  by  adjacent 
grains.  Thus,  they  will  not  exhibit  the  discontinuous  motion  of  the  sand  grain  seen  whoi 
asperities  cm  the  grains  break.  Furthermore,  because  of  their  size,  the  lead  particles  are 
more  uniformly  constrained  than  the  sand  grains.  In  summary,  the  embedded  lead  particles 
provide  insight  into  the  internal  microstructural  response  of  the  specimen  during  loading. 
Limitations  imposed  by  the  use  of  the  particles  tend  to  make  their  response  align  itself  more 
closely  with  displacements  as  predicted  by  the  use  of  a  local  theorem.  The  lead  particles 
may  be  inc2qHd>le  of  experiencing  the  variation  in  mobility  as  is  shown  by  direct  observation 
of  the  sand  grains. 
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Hoarding  the  information  gaiiwd  from  the  epoxy  imi»i^nation  and  sectioning  of  the  tested 
iqsecimen,  general  insight  is  gained  into  the  distributicm  of  damage  (individual  grain  fra^ure 
and  disint^ratiim)  diroughout  the  specimm.  A  qualitative  observation  will  be  presented 
inior  to  riisciiMing  the  types  and  distribution  of  internal  damage.  Damage  is  readily 
apparent  in  the  exposed  specimen  cross  secticms  due  to  a  change  in  the  opadQr  of  individual 
grains.  The  grains  become  opaque  if  they  fracture.  Consider  the  specimens  aibjected  to  the 
Implication  of  10.0-percent  macroscopic  strain.  Damage  is  oxicentrated  at  the  comers  of 
the  specimen  almost  to  the  exclusion  of  damage  lU  the  specimen  int^or.  If  a  sfrimie  of 
some  unit  diameter  encased  within  a  cube  having  sides  of  the  same  unit  length  is  pictured, 
the  portion  of  the  specimen  contained  within  the  sphere  would  be  reladvdy  undamaged,  and 
the  portion  outside  the  sphere  boundary  would  be  relatively  highly  damaged.  Specimens 
which  had  been  subjected  to  various  levels  of  macroscrpic  uniaxial  strain  were  impr^nated 
and  sectioned.  The  damage  appears  to  commoice  at  the  comers  and  propagate  inward  as 
the  magnitude  of  macroscopic  strain  increases. 

Hiese  observations  are  corroborated  by  conclusions  drawn  on  the  displacement  of  embedded 
particles  (Fig.  41)  as  well  as  the  grain  tracks  scored  on  the  ivalls  of  the  confinement  vessel 
(Fig.  40).  It  has  already  been  shown  that  grains  which  reside  adjacent  to  the  vessel  wall  are 
relatively  immobile.  The  immobility  causes  stress  concmtration  on  individual  grains 
resulting  in  fracture  and  disintegration.  This  tendency  appears  to  be  more  prcmounced  in  the 
grains  which  reade  in  the  specimen  comers  (outside  the  boundary  of  the  included  ^here). 
Furthermore,  the  mobility  of  a  given  grain  is  more  highly  restricted  if  the  grain  is  located 
relatively  deeply  within  a  given  comer.  As  the  highly  immobile  grains  located  at  the 
comers  are  cmshed  during  the  initial  stages  of  loading,  the  grain  remnants  are  shifted  into  a 
more  compact  arrangement  and,  thus,  permit  inner  grains  (diose  away  fiom  the  comers)  to 
move  outward  in  the  direction  of  the  comers.  This  activity  occurs  throughout  the  loading 
phase  and,  therefore,  grains  along  the  outer  walls  of  the  confinement  vessel  migrate  toward 
the  comers. 

Grains  located  adjacent  to  the  base  of  the  confinement  vessel,  as  wdl  as  those  in  contact 
with  the  loading  cjq),  show  no  tendency  to  displace  horizontally  during  loading.  This  is 
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readily  explained  by  the  stress  state  in  the  specimen.  In  the  case  of  uniaxial  strain,  a 
vertical  force  is  exerted  on  the  q)ecimen  which  forces  the  grains  to  primarily  displace 
vertically.  Hiere  is  no  external  macroscopic  force  cm  the  grains  causing  them  to  displace 
horizontally.  Grains  in  contact  with  either  the  base  plate  or  loading  cs^  will  generally  not 
displace  horizontally.  Grains  residing  against  the  sidewalls  are  forced  to  displace  vertically 
and  only  exert  a  macroscopic  stress  on  the  sidewall  30  to  40  percent  of  the  vertical  stress. 
Thus,  the  grains  are  relatively  free  to  move  laterally  during  the  ai^lication  of  vertical 
uniaxial  strain.  This  freedom  is  manifested  as  grain  tracks  which  tend  toward  the  comers  of 
the  specimen. 

Figure  42  shows  the  types  of  damage  which  were  ai^)amit  following  epoxy  impregnaticm 
and  sectioning.  The  figure  includes  the  binary  (two-tone)  images  of  two  viewed  sectimis. 
Figure  42  also  shows  other  images  of  damaged  grains.  The  damage  was  generally  limited  to 
fracture  of  grains  with  two  possible  outcomes.  Hther  the  renmants  remained  in  direct 
contact  with  one  another  and  the  grain  appeared  to  still  be  intact  or  the  remnants  broke  zpait 
leaving  a  slight  void  where  the  grain  originally  existed.  Damage  aside,  very  few  of  the 
grains  appear  to  be  in  contact  with  one  another.  This  angle  observation  highlights  the  three- 
dimensional  nature  of  the  connectivity  which  exists  between  discrete  grains.  The  contacts 
exist  either  above  or  below  the  level  at  which  the  sectitm  was  takoi.  Likewise,  damage  has 
a  three-dimensional  nature  which  one  can  only  hope  to  catch  a  glimpse  of  by  employing  the 
technique  of  impregnation  and  sectioning.  It  is  possible  to  come  up  vath  some  measure  of 
the  distribution  of  damage  on  some  arbitrarily  chosen  plane  in  the  q)ecimra.  Likewise,  one 
can  even  quantify  the  type(s)  of  damage.  However,  such  information  can  be  misleading  if 
one  is  trying  to  quantify  the  structure  and  distribution  of  the  damage  in  three  dimensions. 

The  limitation  of  being  able  to  view  the  specimen  (or  section)  in  only  two  dimensimis  is 
again  imposed. 

A  significant  amount  of  damage  became  evident  only  upon  unloading.  Consider  the 
q)ecimens  subjected  to  the  monotonic  application  of  10.0-percent  macroscopic  strain  (Figs. 
43  and  44).  As  mentioned  previously,  individual  grain  fracture  and  disintegration  were 
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(a)  Photomicrograph. 


(b) 


The  OASYS-enhanced  image.  Bonding  agmit  idiown  in  grey. 


Figure  42.  Virgin  ^Mcimoi  (unstrained)  (x  25). 
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(a)  Photomicrogr^h. 


(b)  The  OASYS-enhanced  image.  Bonding  agent  not  shown. 
Figure  43.  Specimen  after  lO-percent  jq^lied  strain  (x  25). 


93 


(b)  Horizontal  cut,  10-percent  macroscopic  strain  (x  25). 

Figure  44.  Typical  images  of  microstructural  damage  (10  percent  macroscopic  strain). 
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(c)  Vertical  cut,  lO-percott  macroscopic  strain  (x  SO). 


(d)  Horizontal  cut,  lO-percent  macroscopic  strain  (x  SO). 
Figure  44.  Concluded. 
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highly  evident  in  the  cinema  records  as  well  as  damage  chains  and  grain  clusters.  However, 
it  ^ipeated  (at  least  quantitatively)  that  an  equal  number  of  grains  disintegrated  during 
unloading  as  were  destroyed  during  initial  loading.  Thus,  fracture  was  much  more  prevalent 
in  the  grains  than  was  2q^>arent  during  loading.  It  was  only  during  the  unloading  jrfiase  that 
these  fractured  grains  made  their  presence  known  by  disintegrating.  This  observation  has 
significant  importance  when  one  considers  the  ^rength  mobilized  by  ^lecimens  subjected  to 
cyclic  loading.  Grain  fracture  occurring  during  the  loading  jrfiasefs)  may  manifest  itself  in 
grain  disint^ration  during  unloading.  The  remnants  then  reorient  themselves  into  mcne 
compact  arrangemoits  during  subsequent  loading.  This  usually  leads  to  the  specimens 
subjected  to  cyclic  loading  exhibiting  lower  md}ilized  peak  strength  at  peak  strain,  as  was 
observed. 
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7.0  CONCLUSIONS  AND  RECOMMENDATIONS 
7.1  CONCLUSIONS 

The  micromotions  of  individual  grains  in  a  q)ecimen  under  uniaxial  macroscq>ic  strain  are 
discontinuous  and  stochastic.  The  magnitudes  of  individual  grain  displacemrats  are  broadly 
distributed  about  some  mean  value  which  is  the  displacement  of  the  "average  grain."  Tliis 
distribution  of  motions  is  exponential  with  relatively  few  grains  being  immobile  within  a 
moving  iwighbortiood.  Some  grains  displace  three  times  as  far  as  the  "average  grain"  for  a 
given  increment  of  applied  macroscopic  strain.  Although  the  underlying  physics  of  this 
micromotion  are  still  under  investigation,  the  nature  of  the  displacement  calls  for  the  use  of 
nmilocal  theories  to  accurately  describe  the  constitutive  resptmse  of  granular  media. 

The  observed  differential  motions  and  evolution  of  damage  may  provide  insight  r^arding 
the  mechanism  of  energy  dissipation  within  granular  material.  Further  insight  will  be 
provided  by  data  obtained  using  an  automated  optical  analysis  system  which  is  capable  of 
providing  highly  consistent  analysis  of  grain  fracture  and  crushing.  The  influence  of  the 
confinement  vessel  walls  (and  ssQ>phire  windows)  on  damage  evolution  requires  further 
investigation  to  correlate  observed  surficial  damage  with  tiie  intoior  damage.  Hiis  is 
because  the  nucrostructural  response  at  the  specimen-container  interface  is  quantitatively  and 
qualitatively  different  from  that  of  the  specimen  int^or. 

Regarding  the  motion  history  of  the  "average  grain"  which  rq)resents  a  neighborhood, 
nonlinear  displacement  fields  have  been  observed  through  direct  observation  of 
micromotions  via  sapphire  windows  and  through  the  use  of  x  rays.  The  displacement  fields 
are  not  directly  proportional  to  applied  macroscopic  strain  as  predicted  by  local  theories. 
Rather,  strain  in  the  immediate  vicinity  of  the  loading  aq>  is  greater  than  the  sq^lied 
macroscopic  strain  with  other  regions  of  the  specimen  generally  experiencing  strain  levds 
less  than  the  applied  macroscopic  strain.  Three  "zones"  masting  within  the  specimoi  under 
uniaxial  strain  are  (1)  a  zone  in  which  relatively  high  local  displacemoits  and  strain  occur, 
(2)  a  zone  bdiind  this  "front"  wherein  the  grains  move  directly  in  relation  to  the  advance  of 
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the  loading  cap,  and  (3)  a  zone  which  has  not  yet  felt  the  influence  of  the  advancing  loading 
csq>.  Sidewall  frictiim,  although  ncM  negligible,  appears  not  to  play  an  overriding  rde  in  the 
grain  micromc^ons. 

A  qualitative  review  of  the  miciomotions  reveals  that  the  grains  appear  to  cmly  shift  position 
and  rearrange  themselves  into  a  more  compact  arrangement  during  the  sqiplication  of  up  to 
4.0-percent  macroscq>ic  strain.  Grain  fracture  and  crushing  become  mote  ai^nrent  at 
higher  macroscopic  strain.  This  process  repeats  itself  with  the  remnants  of  the  fractured 
grains  arranging  themselves  into  compact  orientation  prior  to  their  being  fractured  also. 
These  effects  do  not  randomly  occur.  Each  resprmse  is  observed  in  localized  zones  which 
grow  in  recognizable  patterns.  Fracture  may  or  may  not  occur  before  an  individual  grain  is 
crushed.  Once  fracture  and  crushing  have  begun,  clusters  of  intact  grains  are  perceived  with 
localized  crushing  zones  and  fixture  "chains”  are  evident. 

A  theory  has  been  proposed  which  is  based  on  the  diffusion  equatirm.  This  theory  evolved 
from  the  finding  that  displacement  within  a  neighborhood  of  grains  is  mmuniform  and  obeys 
an  exprmential  distribution.  The  mean  of  the  distribution  is  centered  about  the  motion  of  the 
"average  grain*  which  represents  the  neighbortKX>d.  The  thermodynamic  foundation  of  the 
theory  permits  the  later  incorporation  of  two  properti^  of  granular  media,  configuratkmal 
temperature  and  configurational  entropy.  Both  of  these  prq)erties  are  thought  to  play  a 
central  role  in  the  explanation  of  localization  during  the  deformation  of  granular  media.  The 
material  cmistants  in  this  theory  can,  fot  the  most  part,  be  determined  by  {rfienomenological 
experiments.  However,  microscopic  (microstructuial)  observations  in  three  dimensicms  are 
still  needed  to  shed  light  on  the  underlying  phyncs  of  the  macroscopic  mechanical  behavior 
of  the  media. 

A  brief  description  of  several  classes  of  constitutive  models  for  granular  materials  was 
presented  in  this  report.  The  models  can  generally  be  grouped  undo*  two  headings, 
macroscopic  (phenomenological)  and  microstructural  (particulate).  Both  types  of  models 
have  bem  used  extensively  to  understand  and  predict  the  mechanical  response  of  soil.  But 
each  type  is  limited  in  its  ability  to  predict  the  mechanical  behavior  of  most  naturally 
occurring  soils.  The  macroscopic  n.odels  are  usually  limited  to  predicting  the  response  of 
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soils  subjected  to  well-defined  loading  regimes  within  which  the  general  re^nse  of  the  soil 
is  already  known.  These  models  typically  relate  the  soil  response  to  macroscopic  material 
pr(^)erties  and  do  not  incorporate  the  underlying  physics  associated  with  the  interaction  of 
individual  sand  grains.  Microstructural  models  have  been  developed  to  incorporate  the 
underlying  physics;  however,  these  models  appear  to  be  limited  by  the  current  inability  to 
model  a  sufficient  number  of  irregularly  shaped  grains  of  soil.  This  limitation  is  driv«i  by 
the  need  to  conduct  "experiments"  using  disks  (2D)  and  spheres  (3D)  for  computational 
efficiency.  Experimental  programs  have  been  undertaken  to  explore  the  microstructural 
attributes  of  granular  media  and  understand  how  these  attributes  manifest  themselves  in  the 
material  macroscopic  response.  These  programs  usually  rely  on  viewing  sections  of  soil 
which  have  been  subjected  to  load  and  subsequently  frozen  or  impregnated  with  epoxy. 

These  studies  have  provided  insight  into  what  physical  changes  typically  occur  in  stressed 
soils,  but  they  provide  little  insight  into  the  evolution  of  such  changes.  A  unifying 
technique  is  needed,  therefore,  wherein  one  can  quantify  and  even  view  the  evolution  of 
microstructural  changes  which  occur  in  granular  media  under  arbitrary  loading.  A 
technique  is  proposed  next  which  shows  reasonable  promise  for  relating  mechanical  resptmse 
to  experimentally  observed  microstructural  changes. 

It  is  proposed  that  assemblages  of  irregularly  shaped  particles  be  subjected  to  arbitrary 
loading  using  distinct  element  computer  programs.  Rather  than  using  large  numbers  of 
spherical  or  ovoid  particles,  the  "specimen"  would  be  created  from  the  volumes  rendered  by 
suitable  manipulation  of  images  obtained  from  computed  tomogrtq>hy  (CT)  scanning.  The 
significant  obstacle  would  be  the  creation  of  object-oriented  volumes  from  the  two- 
dimensional  images  provided  by  the  CT  scans.  Determination  of  the  basic  descriptors  of 
fabric  and  its  evolution  under  loading  would  be  trivial.  Furthermore,  advanced  visualizatitm 
techniques  would  permit  one  to  highlight  and  view  only  those  particles  which  undergo 
micromotions  exceeding  preset  thresholds.  Thus,  one  could  isolate  and  study  only  those 
particles  which  participate  in  the  development  of  a  localization  and  which  undergo  excessive 
displacement,  rotation,  or  damage,  etc.  The  specimen  could  also  be  examined  at  any  stage 
during  the  loading  process  without  altering  the  structure  of  the  assemblage.  A  preliminary 
stq)  in  this  direction  has  been  taken  by  integrating  the  use  of  machine  vision  and  artificial 
intelligence  in  quantifying  the  micromotion  and  rotation  of  surficial  particles  of  granular 
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specimens  under  uniaxial  strain.  The  experimentally  determined  microstructural  re^nse(s) 
of  the  surficial  particles  will  serve  as  a  benchmark  against  which  the  microstructural 
reqxmseCs)  of  the  surficial  particles  in  the  distinct  element  computation  can  be  compared.  If 
the  response(s)  match,  it  is  thought  that  the  microstructural  changes  occurring  in  the 
granular  specimen  will  be  mimicked  by  the  changes  occurring  in  the  distinct  element 
specimen.  Thus,  experimental  and  computational  reservations  can  be  linked  to  provide  a 
sound  basis  for  the  developmoit  of  a  microstructure-based  continuum  model  for  granular 
media. 

7.2  RECOMMENDATIONS 

The  difficulty  of  obtaining  consistoit  results  using  manual  methods  during  observatirm  of 
grain  rotations,  fracture,  and  disintegration  dictates  the  need  to  automate  the  data  collection 
process.  This  is  true  in  light  of  the  attempt  to  measure  microstructural  effects  from  the  two- 
dimensional  cinema  images  and  became  more  ^jparent  as  three-dimensional  microstructural 
effects  were  considered.  Suitable  off-the-shelf  image  analysis  programs  are  now  available 
which  can  be  user-modified  to  suit  specific  needs  (filters,  shagpe  fiinctions,  etc.).  Automated 
optical  techniques  should  be  aggressively  explored  (mainly  in  the  area  of  software 
development)  if  three-dimmsional  visualization  and  discrete  element  modeling  of  assemblies 
of  irregularly  shaped  grains  are  to  become  a  reality. 

This  three-dimensional  visualization  and  discrete  element  modeling  should  also  be 
intensively  investigated  as  the  coupling  of  these  techniques  (automated  imaging  and  three- 
dimensional  visualization)  will  aid  in  studying  the  evolution  of  f£d>ric  during  specimen 
loading.  The  specimen  in  this  example  will  be  an  assemblage  of  irr^ularly  shaped  particles 
which  are  in  a  suitable  format  for  manipulation  by  a  discrete  element  program.  The 
primary  activity  in  this  pursuit  will  be  to  investigate  the  existence  of  anisotropy  in  the  initial 
fabric  and  the  evolution  of  both  (anisotn^y  and  fabric)  during  loading.  An  equally 
impcntant  goal  will  be  to  develop  a  deeper  understanding  regarding  the  relationship  between 
these  factors  and  the  occurrence  of  localizations.  These  aspects  are  critically  important  if 
one  is  to  develop  a  three-dimoisional  continuum  model  for  granular  materials  founded  upon 
microstructural  effects. 
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